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PREFACE 

This  report  contains  a major  portion  of  the  lecture  notes  for  a 
graduate  course  in  special  topics  on  advanced  antenna  theory,  that  was 
taught  by  Professor  Giorgio  Franceschetti  of  the  University  of  Naples, 

Italy,  during  his  visit  at  the  University  of  Illinois  at  Chicago  Circle 
in  the  Spring  of  1976. 

The  three  chapters  of  the  report  are  respectively  concerned  with 
antenna  admittance  and  gap  problem,  with  transient  radiation  from 
antennas,  and  with  transform  properties  of  radiated  fields.  Strong 
emphasis  is  placed  on  unconventional  techniques.  Much  of  the  material 
represents  original  research  by  the  Author,  and  is  not  yet  available  in 
textbooks  or  scientific  journals. 

The  members  of  the  Communications  Laboratory  are  grateful  to 
Dr.  Franceschetti  for  his  lectures  as  well  as  for  many  illuminating 
discussions,  and  to  the  Air  Force  Office  of  Scientific  Research  for 
supporting  the  preparation  of  this  report. 


July  1976 


Piergiorgio  L.E.  Uslenghi 
Director  of  the  Laboratory 
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CHAPTER  I 

ANTENNA  ADMITTANCE  AND  GAP  PROBLEM 

"Wherein  the  reader  is  led  down 
the  highways  of  series  manipu- 
lation" 

1 . The  problem 

An  important  problem  in  antenna  theory  in  the  RF  range  is  the  evaluation 
of  the  input  admittance,  i.e.,  the  ratio  of  the  current  to  the  voltage,  at  the 
input  terminals  of  the  antenna  itself.  Motivation  of  this  problem  is  not  only 
academic  in  nature  - the  solution  of  an  elegant  boundary  value  problem  - but 
has  also  deep  roots  in  applications. 

Radiation  of  broadband  signals  requires  a careful  matching  between  the 
output  admittance  of  the  transmitter  and  the  input  admittance  of  the  antenna. 

For  subsurface  applications  the  latter  is  also  depending  on  the  dispersion 
characteristics  of  the  external  medium.  In  the  non-communication  applications 
of  electromagnetics,  measurement  of  the  input  admittance  of  the  antenna,  used 
as  a probe,  allows  possible  determination  of  electromagnetic  constants  of  the 
external  medium,  e.g.,  in  plasma  applications  and  biological  diagnosis. 

The  previous  discussion  highlights  the  importance  of  a rigorous  theory  of 
antenna  input  admittance.  This  is  required  to  give  the  proper  dependence  of 
input  admittance  on  antenna  geometry,  frequency  of  applied  signal,  medium  elec- 
tromagnetic parameters  and  feeding  system  characteristics. 

The  last  one  is  probably  the  most  difficult  point  to  analyze.  Let  us  con- 
sider, in  the  most  general  case,  the  antenna  as  composed  by  two  metal  bodies 
fed  by  a transmission  line  (fig.  1).  The  input  admittance  at  terminals  AA' 
will  depend  on  the  position  of  points  (or  small  areas)  AA'  on  antenna  body  as  “ 

-3- 

L — ■ — — 


r 


well  as  on  the  configuration  of  feeder  structure  close  to  the  antenna  itself. 
Accordingly,  the  input  admittance  cannot  be  specified  in  terms  of  the  antenna 
only.  This  difficulty  can  be  alleviated  by  considering  antennas  with  a sym- 
metry axis  and  fed  at  a small  gap.  This  leads  to  the  idea  of  a 6-gap,  i.e., 
an  infinitesimally  thin  gap  across  which  a finite  voltage  (and,  consequently, 
an  infinite  electric  field)  is  applied.  However,  this  assumption  implies  an 
infinite  value  of  the  imaginary  part  of  the  input  admittance  of  the  antenna. 

It  is  obvious  that  the  above  difficulties  exist  also  in  the  zero -frequency 
limit,  i.e.,  in  the  definition  of  the  input  capacitance  of  the  metal  structure 
of  fig,  1.  This  suggests  the  possibility  of  exploring  the  gap  problem  under 
quasi-static  conditions,  where  the  analytical  difficulties  connected  to  the 
solution  are  smaller.  This  procedure  is  presented  in  the  following  Sections, 
first  with  reference  t the  spherical  symmetry  - a simpler  case  - and  then  to 
the  spheroidal  one.  In  the  limit  of  a very  thin  prolate  spheroid,  the  perfor- 
mance of  thin  wire  cylindrical  antenna  is  recovered. 


Fig.  1 - Antenna  geometry 

The  analysis  which  will  be  presented  leads  to  an  analytical  expression  for 
the  input  admittance  which  clearly  exhibits  the  peculiar  resonant  performance 


M 
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of  antennas.  As  an  application,  transient  input  responses  are  computed  for 
antennas  in  non-dispersive  as  well  as  dispersive  media. 


2.  Relevant  formulas  and  pertinent  expansions 


For  the  readers'  convenience,  relevant  formulas  and  pertinent  expansions 
used  throughout  this  chapter  are  herein  collected. 

Spherical  Hankel  functions  of  the  second  kind 
(note:  the  superscript  (2)  has  been  omitted  for  notational  simplicity) 

h2n-l«  ’^"zn-Z  - h2„-2W 
where  a prime  means  derivative  with  respect  to  z. 


, , , . 2n-2  r . -I  p (2n-2+m)  ! 1 

z h.  _(z)  = l exp[-iz]  ) — ttt — o , , 

2n-2^  ' J £m  m!(2n-2-m)!  ...  .m  * 

0 ' (2iz) 


(2.2) 


where 


^ H„« 


For  m = (z/2)  ' >>1  and  t = [z-  (2n-3/2)  ]/m  > 0 and  not  large 


t H f 2 t3/2  ^ 

j-  ln{ z h (z))  = - - exp(iy)  2/3  ■ | 

dz  2n-2^  m 3 ^^(f  t3/2) 


(2.3) 


For  t-*-0 


2 

E*"Uh2n.2U)]  =■  -;«P(il)(3)1/3^- 

1 


(0.46+i  0.80)  (2.4) 


For  z large  and  z > > 2n-l 


d7  lntZ  h2n-2^  = - 1 


(2.5) 


For  2n-l  large  and  2n-l  > > z 


h2n.2“»  ■ ™ ^ ^ 


For  z-»-0 


2n-l  h2n-l(z) 
2 h2n-2(z) 


2n-l  h2n-l(z) 
2 h2n-2^z^ 


(2n-l)(2n-|) 


1 


Modified  3e3sel  functions  of  the  first  kind 


For  t < < 1 


t n l 

C2  J nT  ’ n - 1 


For  t > > 1 


I (t)  * ®5E££I 


/lit" 


Legendre  polynomials  of  the  first  kind 
Let  n = cos9,  a = tt/2-0 


/ pLi‘'»  pLl"»  d„  . Mitii  sr 


l-"J  pLiw 


P21-1(0> 


2 n n p2n-l(n)  • 2 n P2n(n) 


2 n P2nC0) 


For2n-l>>l,  (2n-l)0>l 


P2n-1^  * 


tt  sin0 


[(2n-i)e  ,J]. 


4n  1 

Tt  ^cos  nlT  cos[(2n-  y )a] 


- <->"  A 


Series  summation 


Accordingly 


* , . , exp(ix) 

I expCinx)  = 

ln 

f dx  l exp(inx)  * l f dx  exp(inx)  ■ j f exP^»£). 
J in  i"  J ln 

•I d*  i!Sp(i»)  • - 1 «"[i-«p(wi 
00 

£ exp(ix)  _ _ £„[i_exp^ix^ ] = 

n ^ 


= - 2,n(2  sin  4 ) - i 


y cos(2n-l/2)x  _ cqs  x 
n 


X y cos_Jnx  + sin 

7 £•  n 


X V 


sin  2nx 
n 


x „ 1 

3 cos  2 *n  7HT7 


-IT  , X 

+ ( ■=■  -X  ) Sin  y 


Integra  Is 


f*> 


dx  = 2 / 1-x  + 


/dx 
x/T! 


1 -X  - 1 


= 2 /1-x  + In 
X /1-x  + 1 


When  the  upper  limit  of  the  integral  is  +1,  compute  the  last  term 
(2.16)  for  x = Xq  and  then  let  xQ-*-l. 


/ 


xv-1(u-x)u"v  exp(xt)  dx  = r ~~  ulJ  ♦ C'W+1;tu) 


where  <}>  is  the  confluent  hypergeometric  function 
4> Cct; 6;  z)  = exp[z]  <p (6-a; 6;  - z) 


(2.14) 


(2.15) 

(2.16) 
of 

(2.17) 

(2.18) 
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Let  us  consider  the  spherical  antenna  of  fig.  2 to  be  referred  to  a 
spherical  system  of  coordinates  (r,n=cos9,<{>)  centered  at  the  center  of  the 
antenna.  The  antenna  is  assumed  to  be  fed  by  a 4>- independent  voltage  Vq 
exp(iuit)  at  a small  finite  equatorial  gap  of  thickness  2s,  s=aa q.  The  external 
medium  is  assumed  linear,  homogeneous,  isotropic,  time- invariant  and  described 
by  permittivity  e and  permeability  u. 

I 

P(r,9) 


Fig.  2 - Geometry  of  the  spherical  antenna 

The  boundary  value  problem  associated  with  the  radiation  from  the  antenna 
is  the  solution  of  steady-state  source-free  Maxwell  equations  subject  to 

(i)  radiation  condition  at  infinity 

(ii)  symmetry  condition  relative  to  ^-coordinate 

(iii)  boundary  conditions  at  r = a. 

Electromagnetic  fields  specified  in  a spherical  coordinate  system  and 
subject  to  condition  (ii)  decouple  into  two  sets  of  modes  characterised  by 


m 


ft 
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components  (H  ,H 


For  the  assumed  excitation  which  has  no 


^-component  of  the  electric  field,  the  former  modes  are  not  excited.  Accord 


ingly,  the  solution  of  Maxwell's  equations  subject  to  condition  (i)  in  spheri 
cal  coordinates  is  the  following: 


where  h (x)  is  the  spherical  Hankel  function  of  the  second  kind,  P1(x)  the 
m m 

associated  Legendre  polynomial  of  the  first  kind,  k = u> /ey",  and  t are  as  yet 

n 

unknown  expansion  coefficients.  Note  that  the  usual  superscript  (2)  in  the 
notation  of  the  spherical  Hankel  function  has  been  omitted  for  the  sake  of 


notational  simplicity.  Forcing  of  boundary  condition  at  r = a requires 


E being  the  9-component  of  the  applied  electric  field  at  r = a.  Note  that  E = 0 


If  V(r|)  is  the  applied  voltage  distribution,  the  applied  field  distribu 
tion  E(n)  follows  immediately  from  (3.4): 


The  expansion  coefficients  tn  appearing  in  (3. 1-3.3)  can  be  computed  by 


,1 


(i)  evaluating  (3.2)  at  r=a, 

(ii)  multiplying  both  sides  of  (3.2)  by  P^m  ^(n) 

(iii)  integrating  in  n-domain  between  (-1,1) 

(iv)  taking  advantage  of  orthogonality  relation  (2.11). 
Hence: 


icoe 


4n- 


Cn  (2n-l)  h2n  j(ka)-kah2n  2(ka)  4n(2n 


/n0 

/I^n7  v'(n)  p2n.i(n)dn  (3.6) 

'no 


Let  us  define  the  input  admittance  Y as  the  ratio  of  the  current  at  the 
gap  to  the  applied  voltage 


2ira  H^(a,n0) 


Using  (3.3  and  3.6)  into  (3.7)  the  formal  result 


(3.7) 


Y = iiruiea  £ 


4n-l 


R (ka) 


V (nJ 


xn  (2n-l) 2 "h^-1 


(2n-l)  h.  .(ka) 

R (ka)  = - u u - rt  — 

n ka  h2n  2(ka) 


VV  ■ - 


fr^-f  ° p2n-i^J<in 

0 


(3.8) 


(3.9) 


(3.10) 


is  obtained. 

The  following  comments  are  now  in  order. 

(i)  The  expression  V^n^)  can  be  considered  as  the  normalized  excitation 
voltage  of  the  (2n-l)  mode.  For  a 6-gap  excitation,  i.e.. 


V'(n)  = - Vg6(n)»  the  result 
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(3.11) 


(P2n-l(0)l  2 

vn 2n  lp2n(°)) 


is  obtained.  Note  that  use  has  been  made  of  (2.12). 

(ii)  For  large  values  of  n and  ka  fixed,  eqs . (2.7  and  2.14)  show  that 

8n2 

R (ka)  a — j , V - it . Accordingly,  the  series  (3.8)  diverges  for 
z 

a 6-gap  excitation. 

(iii)  Since  this  divergence  process  does  not  depend  on  ka,  it  is  conve- 
nient to  add  and  subtract  to  the  series  (3.8)  its  value  for  ka  = 0. 
Using  (2.8),  the  result 

Y (ka)  -Y(0)  + [Y (ka)  - Y (0) ] = 


3 iumea  £n  W + 

in  (2n-l) 


r 4n-l  1 „ , . 

iwirea  £ j Okal-l  W 

xn  (2n-l)  VKaJ  1 n u 


(3.12) 


is  obtained. 

The  first  series  in  (3.12)  is  still  diverging  for  a 6-gap  excitation. 
However,  the  second  series  is  now  converging  and,  for  nQ  < < 1,  the  6-gap  exci- 
tation is  a reasonably  good  approximation.  Accordingly, 


Y = iwTrea  f — — — =-  V (n0)  + 
l"  (2n-l) ^ n U 


+ iwrrea  f [ P (0)  ]2  r-*T- - 

(2n-l)  2n  V*  J 1 


(3.13) 


Aside  from  the  factor  i,  the  first  term  in  (3.13)  is  clearly  the  susceptance 
associated  with  the  static  capacitance  of  the  antenna.  Its  value  will,  in  gen- 
eral, depend  on  the  gap  dimension  (nQ)  and  type  of  excitation  (E(n)).  The  second 
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series  approaches  zero  as  w-5  for  (see  eq.  (2.8))  in  a lossless  medium. 

Accordingly,  it  is  a typical  dynamical  term  and  will  be  called  the  dynamical 
input  admittance  Y^. 

The  formal  expression  of  the  input  admittance  is  then  the  following 

Y = iw  Cs  + Yd  (3.14) 

4 . The  spherical  antenna  continued.  Static  capacitance 

From  (3.13)  the  expression  of  the  static  capacitance  of  the  antenna  is 
given  by 

Cs  = irea  l i"---1  - V (n  ) (4.1) 

3 !n  (2n-l)2  n 0 

In  order  to  sum  (4.1)  the  distribution  of  the  field  V (n)  applied  across 
the  gap  should  be  known;  this,  in  turn,  requires  specification  of  the  feeder 
geometry  and  the  solution  of  a generally  complicated  boundary  value  problem. 

It  is,  therefore,  customary  to  assume  a priori  the  voltage  distribution  across 
the  gap;  this  is  acceptable  only  if  Cg  is  rather  insensitive  to  the  field  dis- 
tribution across  the  gap.  This  is,  however,  a rather  crucial  point  on  which 
the  very  definition  of  an  input  admittance  of  the  considered  antenna  must  be 
founded.  As  a matter  of  fact,  this  definition  is  possible  only  on  the  assump- 
tion that  Y is  (practically)  independent  of  the  field  distribution  across  the 
gap  and  therefore  does  not  (practically)  need  the  feeder  specification. 

Let  us  define 

oo  4n  - 1 oo 

W = -772  VV’  <4-2> 

(2n-l) 

00 

where  is  the  asymptotic  expansion  of  where  n is  very  large.  The  expres- 
00 

sion  for  f can  easily  be  computed  by  using  the  asymptotic  expansion  (2.13). 
Hence: 


<»  pLi‘V  pLi<V  ■-  r 


{cos  [ (2n  -i)(a+a  )]  + 


cosOq  cosa 


+ cos  [(2n  - j)(a-aQ)]}  . 


- r: 


cos [m(a-aQ) ]da 


cos[m(a+aQ)]da; 


this  result  stems  from  the  change  of  variable  a-*- -a,  and  remains  valid 
when  an  even  function  of  a is  introduced  inside  the  integral. 


Un)  vn(n0) 


TT/COSa 


/an  , 

/cosa  V ’ (a)  cos  [(2n  - j ) (a+a^)  ]da  . 


0 -a 


Accordingly, 


W = 


~2  ■ ■■■  / / cosa  cos[(2n  - j ) (a+a.)]da  . 

iTn  /cosaQ  -cIq  u 


Now  the  series  (4.1)  can  be  recast  in  the  following  form: 


Cg  = nea 


l ^ * *“  [jn  ^ W • *V] 


(4.3) 


Note  that  in  the  second  series,  each  term  behaves  as  2/ffn  plus  terms  of 
2 

order  1/n  as  n-*-°°.  Accordingly,  the  series  is  uniformly  convergent  with 
respect  to  nn,  and  its  evaluation  can  safely  be  performed  for  nn  = 0.  Hence 


l j 2n(4n-l)  f 2 _ 2 _/  al>128 

(2n-l)2  2n  1fn( 


(4.4) 


The  first  series  which  appears  in  (4.3)  is  also  uniformly  convergent  with 
respect  to  nn,  thus  allowing  interchange  of  summation  and  integration  opera- 


tions. Using  (2.15): 
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I f>0>  ■ 


it  / cosa 


0 -«r 


2 sin(a+aQ) 


0*+aJ 


1 • a+an 
J Sln— I 


'cosa  da  = 


/a 


0 V (a) 


2(a+aQ)  * 


(4.5) 


the  last  result  being  valid  for  aQ  <<  1. 

Assume  now  V^-a^)  =0;  V(a)  = -Vq.  Integration  by  parts  of  (4.5)  leads  to 


l f°°(nn)  Z*1  dx 

,n  n 0 * 4a  it  J x 


v(0)  = 0 ; vCl)  = 1 


(4.6) 


Clearly,  only  the  second  term  of  (4.6)  depends  on  the  voltage  distribution 
across  the  gap. 

Summing  up  (4.4  and  4.5),  the  result 

Cc  = 27rea  f 0 . 564  + - In  §-  + - f dxl  = 

s L IT  4S  TT  J xj 

0 

= 2irea  [0.882  + ^ JJ  + I (v)  ] = 2ea  [Jin  + I(v)],  (4.7) 


with  I(v) 


■fl 


1 dx 


(4.8) 


is  obtained.  Clearly,  only  the  term  I(v)  depends  on  the  voltage  distribution 
across  the  gap. 

An  idea  of  the  magnitude  of  the  quantity  I(v)  can  be  obtained  in  two  limit- 
ing cases.  For  a constant  applied  electric  field  distribution  across  the  gap 
v(x)  =x  and  I(v)  =0.  For  a gap  with  sharp  rims  the  electric  field  should  show 
a square-root  singular  behaviour  at  the  gap  edges.  Accordingly,  a reasonable 
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1 


normalized  voltage  distribution  is 


v(x)  = 


_ 

2 


0 < x < | 


v(x)  = 1 - P 5 , 4 < X < 1 


Using  (2.16)  for  evaluating  the  integral  (4.8)  the  result 
I(v)  = 1.404 


(4.9) 


is  obtained.  Accordingly,  the  static  capacitance  will  be  practically  insensi- 
tive to  applied  field  distribution  as  long  as 


In  — > > 1.404 
s 


(4.10) 


For  a/s  = 10,  i.e.,  ag  = 5.75°,  the  change  in  the  static  capacitance  in  the 
two  considered  cases  is  as  high  as  32%,  for  a/s  = 100,  i.e.,  ag  = 0.57°,  it  is 
still  21%.  The  change  becomes  negligible  for  the  case  of  wire  antennas  (see 
Sec.  7). 

5.  The  spherical  antenna  continued.  First  term  of  the  dynamical 
input  admittance. 

From  (3.13)  the  first  term  of  the  series  representing  the  dynamical  part 
of  the  input  admittance  Y is  given  by 


1 


Yl  = ia)£a  2 R.  (ka)-l 


hj(ka) 

Rl(-ka)  = ka  hQ(ka) 


1 + ika 
(ka)2 


(5.1) 

(5.2) 


where  h^  and  h^  can  be  computed  via  (2.1  and  2.2). 
Substituting  (5.2)  into  (5.1)  the  result 


v 3rr  . 

Yj  ■ y~  1U)£a 


(ka)4  3tt 
' = 2 


1+ika- (ka) ' 


(uea) ‘ 


1 


iuya 


— + ltnea 


3 tt  , .2 

— (we  a) 

/e?U  yx 


(5.3) 
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MM* 


is  obtained. 


Expression  (5.3)  has  the  lumped  equivalent  circuit  representation  given 
in  fig.  3,  where  the  box  represents  a reciprocal  gyrator  such  that  the  input 
admittance  Y^  equals  Yq  divided  by  the  output  admittance  / e/y  y^ . This 
gyrator  is  the  dual  of  a conventional  transformer  wherein  the  mutual  inductance 
is  replaced  by  a mutual  capacitance  / 3 it/ 2 ea  and  terminating  inpedances  by 
terminating  admittances. 


L*ya 


Fig.  3 - Equivalent  circuit  for  the  first  term 
of  the  dynamical  input  admittance. 


The  circuit  of  fig.  3 correctly  represents  the  first  resonance  properties 
of  the  antenna.  The  normalized  admittance  y is  exactly  a parallel  resonant 
circuit.  Due  to  the  presence  of  the  gyrator,  the  network  behavior  of  the  first 
resonant  mode  of  the  antenna  is  a sort  of  inverted  resonance. 

The  admittance  Y^  is  purely  real  when  y^  is  purely  real,  i.e.,  for 
u^/iya  » 1 in  a lossless  medium.  At  this  angular  frequency,  which  is  identified 
as  the  first  dynamical  resonant  frequency  of  the  antenna. 


Y 


1 


For  oj  < < 


Y 


1 


3 


3tt 

2 /y 


(5.4) 


the  term  1/iujya  is  prevaling  in  the  denominator  of  (5.3)  and 


3 

z 


(S.5) 
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where  z=oji/Eua  can  be  considered  as  a normalized  frequency.  The  contribution 
of  this  mode  to  the  radiation  process  becomes  negligible  and  the  only  effect 
is  the  lowering  of  the  (capacitive)  static  reactance  iu)Cs. 

For  u>  > > the  prevaling  terms  is,  on  the  contrary,  icoea  and 

Y1  3 - 1 2 /jT  z t5,6-* 

Again  the  contribution  of  this  mode  to  the  radiation  process  is  negligible  and 
only  reactive  local  fields  are  excited.  These  results  are  valid  also  for 
higher  order  modes,  as  shown  in  Sec.  6. 

In  the  case  of  a dispersive  medium,  i.e.,  e = e(u>)  and/or  u=y(c*0,  the 
above  conclusions  must  be  modified.  However,  the  theory  just  presented  allows 
for  an  easy  study  of  the  input  properties  of  the  antenna  in  such  media.  For 
instance,  in  the  case  of  an  imperfect  conductor  of  conductivity  a,  it  is  only 
necessary  to  substitute  e+a/iw  to  e into  (4.7)  and  (5.3).  This  is  a great 
advantage  with  respect  to  other  antenna  input  admittance  theories. 


The  spherical  antenna  continued.  Higher  order  terms  of  the  dynamical 
input  admittance 

The  general  term 


Y"  ■ iuT£a  ^ [P2"C0”2  ^ 


of  the  series  appearing  in  (3.13)  is  now  considered. 

Expression  (6.1)  can  be  recast  in  the  following  form 


2n(4n-l)  _r„  2_2  2 

( 2n'-~ L ) * fftP2n(0)]  * e a 


(2n-l)  h (z) 

y = iZ  + r-r 

1 h2n-2(z) 
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Now  (2.1)  is  used  to  transform  (6.3)  into 

(2n-l) 2 . ..  d 

y * 12  ♦ - — : — — ♦ i(2n-l)  — 


1 C2n~ 1 ) In  [z  h2n_2(z)] 


Eq.  (6.4)  shows  that  /e/u  y^  can  be  represented  as  a parallel  connection 

2 

of  three  terms:  a capacitance  ea,  an  inductance  ya/(2n-l)  and  a residual 

admittance  /e/y  yf  given  by  the  last  term  of  (6.4),  as  shown  graphically  in 
fig.  4. 


iFya- 


/Fyr 


Fig.  4 - Lumped -circuit  representation 
of  the  term  y . 

Use  of  (2.2)  simply  shows  that  the  normalized  residual  admittance  yr  is 
the  parallel  connection  of  a normalized  conductance  (2n-l)  and  a normalized 
admittance  given  by  the  ratio  of  two  polynomials  of  degrees  (2n-2)  and  (2n-l) 
respectively. 

Clearly,  the  equivalent  lumped-circuit  representation  of  (6.2)  is  a 
reciprocal  gyrator  closed  on  the  normalized  admittance  (6.4).  Accordingly, 
the  equivalent  circuit  for  the  total  input  admittance  Y of  the  antenna  is 
given  in  fig.  5,  where 


,2  2n(4n-l) 


2 2 2 2 


(2n-l) 


j *[P2n(0)r  oTe‘a 


The  above  results  are  exact.  However,  in  order  to  get  a physical  insight 
into  the  resonant  properties  of  the  antenna,  it  is  convenient  to  have  approxi- 
mate expressions  for  y . This,  in  turn,  corresponds  to  an  approximate  evalua- 
tion of  the  last  term  of  (6.4)  in  different  ranges  of  z. 
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Fig.  S - Equivalent  network  for  the  input 
admittance  Y of  the  antenna. 


We  will  restrict  our  attention  to  large  values  of  z.  For  z = 2n-l,  eqs. 
(2.3  and  2.4)  show  that 


yr  3 


2n-l 


1Z 


Yjj  (0.46  + i 0.80) 


Accordingly 


fE  _ , A on  fz  2n-l  (2n-l) 2 f\  . 0.46  1 

/ W n / y T7T  luua  (2n-l)z1/3J 

' r,  0.46  1 

r 


= icoea  ♦ 0.80  (2n-l)2/3  ♦ C2n--1)  | 1 ♦ 

/ y v iwua 


(6.6) 


The  network  behavior  of  the  antenna  is  clearly  that  of  a parallel  resonant 


circuit  which  resonates  at  a frequency  slightly  larger  than  2n-l. 


For  2 <<  2n-l,  taking  (2.6)  into  account 


e 2 ,,  ,.2 

-Yn  XJ  (2n-l)  . 


while  for  z >>  2n-l , taking  (2.5)  into  account 


— y = iz  . 
U 7n 


(6.7) 


(6.8) 


Accordingly,  also  for  the  general  mode  the  network  behavior  of  the  antenna 

is  that  of  an  inverted  resonance.  Each  mode  is  strongly  excited  only  in  the 

neighbor  of  its  resonance  at  the  normalized  angular  frequency 

z-  . = w_  , /iu  a = 2n-l . 

2n-l  2n-l 

7 . The  wire  antenna 

The  results  shown  under  Secs.  3 through  6 are  qualitatively  valid  also  in 
the  case  of  spheroidal  antennas,  although  the  mathematical  machinery  becomes 
much  more  complicated.  However,  the  techniques  of  analysis  already  introduced 
can  be  successfully  applied  also  to  this  case.  Hereafter,  only  results  for 
thin  prolate  spheroidal  antennas  are  quoted,  i.e.,  for  antennas  characterized 
by  large  values  of  the  parameter 


n=2  2,n|i-2  = 2^n  £^36a  (7.1) 

b b 

where  2a  is  the  larger  and  2b  the  smaller  axis  of  the  prolate  spheroid.  The 
obtained  results  can  then  be  applied  to  wire  antennas  of  length  2a  and  diameter 
2b  wherever  (2  > > 1. 

The  equivalent  network  circuit  for  the  input  admittance  of  the  antenna  is 
still  of  the  type  represented  in  fig.  5.  Particularly,  for  the  first  few  terms 
we  have 

Cs  = csea,  (7.2) 
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He*  Is  t*"r‘ I(v)1}  • 


v2  3-rr  2 2 2 

Yoisr  u£a  * 


(7.4) 


'e  fz  , 1 

— y,  = luieac  ♦ / — b + -i r-  , 

U 71  / U lojpaJl  ' 


(7.5) 


where: 


c = 


F 1 + 0.0826  — "1 

L (n  + TT)J 


* n , 


(7.6) 


100  25  1_ 

8 99  11  ,2 

, , 11  ± 

25  n * 


(7.7) 


(7.8) 


, n + 5/n 

* = — — - i • 


(7.9) 


Eq.  (7.3)  shows  that  the  second  term  in  the  curly  brackets  is  negligible 
with  respect  to  the  first  one  when  b/a  is  very  small  ( thin  antennas ) . Under 


this  assumption 


(7.10) 


and  this  static  capacitance  becomes  practically  independent  from  gap  thickness 
and  feeder  characteristics.  This  is  the  very  condition  under  which  an  input 
admittance  can  be  properly  defined  and  also  explains  the  reason  why  theories 
which  do  not  take  into  account  the  gap  problem  give  acceptable  results  for  thin 
wire  antennas. 

In  the  low  frequency  range,  when  only  the  first  mode  of  the  antenna  is 
excited,  the  input  admittance  of  the  antenna  is  given  by 


Y = iooC  + Y = A p-  iz  .9z--iMz  

U ' 9z  + i£ty[(2z)  v-9] 


(7.11) 
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4 


The  antenna  resonates  when  Y is  purely  real.  For  large  values  of  ft  this 
. 2 2 

happens  for  (2z)  \p-9  = 0 and  for  z tp-9  = 0 , as  easily  follows  by  neglecting  the 
term  9z  respectively  in  numerator  and  denominator  of  (7.11). 

In  the  first  case 


z 


1 


3_ 

2<p 


i.e. , a=  0.955  j 


(7.12) 


\ being  the  wavelength:  = 2tt/\ • This  is  the  well-known  first  resonance 

of  the  antenna.  The  input  conductance  at  this  resonance  is  given  by 


i 


G = 


- 

y 2 


(7.13) 


which  correspond  to  the  usual  80  ohm  of  input  resistance  at  resonance. 

The  second  resonance  occurs  at  z = 3,  i.e.,  in  a frequency  range  in  which 
the  contribution  of  the  second  higher  order  mode  is  no  longer  negligible. 
Accordingly,  previous  hypotheses  are  no  longer  valid.  Furthermore,  the  conduc- 
tance  at  this  resonance  turns  out  to  be  proportional  to  l/ft",  so  that  this 
resonance  is  of  no  practical  interest. 


8 . Transient  response  of  the  antenna.  Free-space  environment 

Let  us  now  consider  the  transient  input  current  response  I (t)  of  the  wire 
antenna  when  a time-varying  voltage  V(t)  is  applied  to  the  input  terminals.  It 
is  explicitly  assumed  that  the  frequency  content  of  V(t)  is  negligible  in  the 
high  frequency  range,  so  that  the  input  admittance  Y of  the  antenna  can  be 
approximated  by 

Y a iioCs  + Y:  . (8.1) 


Let  Cq  and  be  the  permittivity  and  permeability  of  free  space,  and  let 
T = /EqPq  a.  If  the  applied  voltage  is  a pulse  of  time  duration  T,  its  frequency 
content  is  of  order  1/T.  On  the  other  hand,  assumption  (8.1)  requires 
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I 


id  /e„un  a = oj  t « 2 . 
max  00  max  ^ 


(8.2) 


Accordingly, 


|2L<1 

T - T 


i . e. , T > ttt 


(8.3) 


Note  that  2t  is  the  time  required  for  the  signal  to  travel  a distance 
equal  to  the  length  of  the  antenna.  Accordingly,  the  pulse  duration  should 
exceed  this  time. 


P = ia)/£nPn  a 


(8.4) 


be  the  analytical  continuation  of  the  normalized  angular  frequency  ou*^e: Q a,  and 


. f 

2-iti  J 


r /£o/Mo 


exp(pt)  dp 


(8.5) 


be  the  normalized  unit  response  (i.e.,  the  current  response  to  a voltage  exci- 
tation 5(t))  of  the  approximate  circuit  corresponding  to  (8.1).  In  (8.5),  the 
integration  path  T runs  parallel  to  the  imaginary  axis  leaving  to  its  left  all 
the  singularities  of  the  integrand,  and  t is  the  time  normalized  to  t. 

The  input  current  response  I(t)  of  the  antenna  follows  by  superposition  as 


l(t)  = /jT  J v(t,)  dt’  * 

0 _co 


(8.6) 


where  normalized  times  are  used. 

In  the  general  case  the  environment  of  the  antenna  is  time-dispersive, 
i.e.,  P=pQpr(p)  and  e=£0er(p),  where  the  relative  permeability  Pr(p)  and  the 
relative  permittivity  er(p)  describe  the  dispersive  properties  of  the  medium. 
It  is  also  convenient  to  introduce  the  relative  refractive  index  nf(p)  = 

= </£r(p)Ur(p)  • 
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By  substituting  the  above  into  (7.11)  and  neglecting  terms  of  order  1/(1 


with  respect  to  unity,  the  result 


is  obtained.  Note  that,  insofar  as  the  third  term  in  (8.7)  is  concerned,  the 


integration  contour  T can  be  closed  by  means  of  a circle  of  infinite  radius  in 


the  left-hand  plane  for  t>0,  so  that  the  integral  is  simply  given  by  poles  and 


branch-cut  contributions 


Let  us  now  consider  the  case  of  a free-space  environment,  i.e 


In  this  case  no  branch-cut  contributions  are  present,  but  only  two  poles  at 


Accordingly,  the  integral  (8.5)  can  be 


s incu  t ] 


Eq.  (8.8)  can  be  recast  in  the  same  approximation  under  the  following  form 


Eq.  (8.9)  clearly  shows  that  the  unit  input  current  response  of  the 


antenna,  for  t > 0,  equals  that  of  a resonant  damped  circuit  with  resonance 


frequency  wn  and  with  damping  constant  an.  Obviously,  the  damping  of  the  cur 


rent  has  a counterpart  in  the  radiation  losses  of  the  antenna 


9.  Transient  response  of  the  antenna  continued.  Dispersive  environment 

In  the  case  of  a dispersive  medium,  the  more  general  expression  (8.7)  has 
to  be  substituted  into  the  integral  (8.5).  Accordingly,  the  following  changes 
are  generated. 

(1)  The  locations  of  the  poles  in  the  complex  p-plane  are  shifted,  this 
taking  into  account  coupling  of  the  antenna  with  the  dispersive 
environment. 

(2)  Branch-cuts  should  be  made,  due  to  the  presence  of  the  square-root 

nr(p)  - /er(p)ur(p)  • 

For  example,  for  a conducting  medium  of  conductivity  a, 

er(p)  « 1 + £ , Ur  - 1 (9.1) 

where  y = aa/Ze^/u^  is  the  normalized  conductivity,  and  now  is  the  dielectric 
constant  of  the  external  medium  (different  from  that  of  free-space) . 

The  branch-cut  is  along  the  negative  ci-axis  from  a = 0 to  a = -y  (see  fig.  6). 


Fig.  6 - Plot  of  poles  and  branch-cuts  in  the  complex 
p-plane  for  the  case  of  a conductive  medium. 
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4 


Let  us  first  consider  what  happens  to  the  pole  contributions  due  to  the 

zeros  of  the  denominator  of  the  third  term  in  (8.7). 

Let,  as  in  the  case  of  free-space,  p^  = -aQ  + iuig,  p2  = -aQ  - ioo0; 

2 

aQ  = 9/8S1  i|T,  u<0  = 3/2/J7.  In  the  case  of  a conductive  environment,  the  poles 
are  possible  solutions  of  the  following  equation 

(^r  P-P^C^j.  P-P2)  = 0,  (9.2) 


where  er  is  given  by  (9.1).  The  solutions  of  (9.2)  are 

P'  = - \ + /[-a0±i(w0  + *)]  [-aQ  ± i(a)Q  - J)]  = -a^  ± iu>'Q  (9.3) 


where  the  plus  sign  has  been  chosen  for  the  square-root  in  order  to  recover  the 
solution  of  Sec.  8 for  y = 0. 

The  pole  contribution  i^(t)  t0  the  total  current  response  (8.5)  can  be 
easily  computed  as 


TT  2 M0 

y = zn  exp(-y)  coswQt  + — ),  (t > o) , 


(9.4) 


which  is  the  generalization  to  the  case  y t 0 of  the  analogous  term  appearing 
in  (8.9).  Accordingly,  is  recognized  as  the  new  resonant  frequency  of  the 
antenna  and  as  the  new  damping  constant. 

The  following  observations  are  now  in  order. 

(i)  For  y < < 2Wq  , eq.  (9.2)  simplifies  to 

P'  = - \ - ± iu>0.  (9.5) 

Accordingly,  only  the  damping  constant  is  increased  with  respect 
to  the  free-space  case,  whereas  the  resonant  frequency  remains 
the  same. 

(ii)  For  increasing  values  of  y the  damping  further  increases  and  the 
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resonance  frequency  diminishes.  There  is  a critical  value  yc  = 2u)q 
of  the  conductivity  for  which 


P’  3 * wf 


'/^Fo 


± 1 


(9.6) 


(iii)  As  the  conductivity  is  further  increased,  the  damping  increases  and 
the  resonant  frequency  of  the  antenna  is  still  lowered.  For 
y > > 2u>q,  we  have 


P' 


-Y 


2°o“o 


(9.7) 


Let  us  now  consider  the  branch-cut  contributions  due  to  the  square-root 
of  £r(p)  which  appears  in  the  first  and  third  term  of  (8.7).  Note  that  the 
closure  of  the  T contour  of  integration  in  the  left-hand  plane  requires  the 
subtraction  from  the  first  term  in  (8.7)  of  its  limiting  value  for  |p|  ■*  *>. 
Then,  letting  x = -a,  the  expression  that  the  sum  of  these  two  terms  assumes 
on  the  upper  side  of  the  cut  is  the  following: 


2 * 

-x(y-x)  -2iaQ/x(Y-x)  + wQ 


(9.8) 


and  is  equal  to  the  complex  conjugate  on  the  lower  side  of  the  cut. 

Summing  up  the  two  contributions  along  the  two  sides  of  the  cut  and  inte- 
grating between  (0,-y)  gives  2iri  times  that  part  of  the  current  response  (8.5) 
which  is  due  to  the  branch  singularity.  Letting  ( t.)  be  such  contribution, 
we  have 


■ £ / 


_ 3/2  “0+*(Y-x) 

Sx  (y-x)  — = y exp(-xt)dx 

[w0  * x(y-x)] 


(9.9) 


where  terms  of  order  (a0/2cj0)^  have  been  neglected. 
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For  coQ  > y/2  we  can  safely  neglect  the  term  x(y-x)  with  respect  to  coQ. 


Then,  the  integral  can  be  computed  by  repeated  use  of  (2.17)  through  (2.20). 
Hence 


3 

iK(t)  » — * Y exp( 


Yt.  r vFO 

- r3  L lt7T  + “Yt72“J 


(9.10) 


For  Yt/2>>1,  eq.  (2.10)  transforms  (9.10)  into 


4/fr 


w. 


$ 


3/2 


(9.11) 


Eq.  (9.11)  shows  that  the  branch-cut  contribution  is  the  dominant  term  at  large 
observation  times.  This  is  a general  result  whenever  the  branch-cut  has  at 
least  one  point  on  the  imaginary  co-axis. 

The  remaining  contributions  to  the  total  unit  current  response  are  straight- 
forward to  compute.  The  final  result  is: 


i(t)  = ~i  fi'(t)  ♦ (6a0+Y)5(t)  ♦ ip(t)  ♦ ib(t). 


(9.12) 


CHAPTER  II 


TRANSIENT  RADIATION  FROM  ANTENNAS 

"It  is  worth  emphasizing  that 
our  knowledge  about  the  tran- 
sient response  of  antennas  is 
very  meager  indeed" 

(T.T.  Wu , in  Antenna  Theory) 


1 . The  problem 

Although  the  steady-state  radiation  from  antennas  has  been  widely  studied 
in  the  past  and  the  radiative  mechanism  is  well  understood,  the  radiation  under 
transient  conditions  is  still  an  almost  undiscovered  field.  In  spite  of  this, 
the  problem  is  rather  important  from  both  theoretical  and  practical  points  of 
view. 

In  order  to  complete  our  knowledge  in  the  field  of  antennas,  it  is  certainly 
highly  desirable  to  have  a complete  understanding  of  radiation  properties  under 
steady  state  and  transient  conditions  as  well.  On  the  other  hand,  wide  band 
transmission  of  signals  and  pulsed  propagation  in  non-dispersive  and  dispersive 
lossless  and  lossy  media  (e.g.,  subsurface  radiation,  radio  links  between  two 
satellites  in  the  ionosphere)  requires  a sharp  understanding  of  the  mechanism 
of  transient  radiation  in  non-dispersive  and,  even  more,  in  dispersive  media. 

It  is  worth  emphasizing,  however,  that  the  analytical  machinery  necessary 
to  correctly  describe  transient  radiation  is  rather  formidable.  Accordingly, 
it  is  very  likely  that  the  physics  of  radiation  could  be  submerged  into  an  ocean 
of  difficult  mathematics.  Therefore,  it  seems  attractive  to  develop  simple 
models  which  could  help  in  the  understanding  of  the  physical  mechanism  of  radia- 
tion with  a minimum  of  mathematics. 
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In  the  following,  first  some  heuristic  procedures  are  presented  which  allow 
for  an  easy  understanding  of  the  radiation  from  elementary  sources,  coaxial  aper- 
tures, wire  and  loop  antennas  in  free-space.  The  problem  of  elementary  sources 
in  dispersive  media  is  touched  upon  as  well.  Then,  a rigorous  solution  of  radia- 
tion from  an  antenna  of  simple  geometry  - the  spherical  antenna  - is  presented. 


2.  Relevant  formulas  and  pertinent  expansions 


Properties  of  the  Dirac  function 


6(at+b)  - i 6Ct  ♦ |) 


(2.1) 


t 5 (t)  = - 6(t)  ; t*  6 (t)  « - 2 t 6 (t) 


(2.2) 


Inverse  Laplace  transforms 

(Fourier  transforms  are  obtained  by  letting  p = icj) 


F(p)  exp(-pt) 


f(t-T) 


I x (px)  exp(-px) 


x - t 


for  -x<t<x,  and  zero  for  j t [ > x 


(2.3) 


IQ(px)  exp(-px) 


TT  J 2 2 

/x  - t 


(2.4) 


for  -x  < t < x,  and  zero  for  |t|  >x 


Ij(PX) 


exp(-px) 


1 y2  72 
— /x  - t 
TTX 


for  -x < t < x,  and  zero  for  J 1 1 > x 


i 


In(px)  exp(-px) 


c-ir  i 


cos 


[n  cos'1  i] 


^ /"3  f 

/x  - t 

for  -x  < t < x,  and  zero  for  |t|  >x 


(2.6) 


In(px)  exp(-px) 
(px)2  ♦ a" 


(-D 

irax 


- C c°-d“  c°s~-  u/x.l.  Sin  r i (t-u)  1 

J.y,  fl  I L X J 

A / x - u 


du  (2.7) 


Bessel  functions 


-X  < t < x 


px 


In(px)  exp(-px) 
(px)2  >a2 


(~i)n  r 

iTax  / 
*'-x 


1 


cos [n  cos  u/x] 

fl  2 
/x  - u 

-X  < t < X 


cos 


[i  (t-U)] 


du  (2.8) 


/~2  7 J , [x /t(t^2r)  ] 2 

exp[-T  /p  + x“  + xp]  6 (t)  - x t U(t) 

x / 1 (t+2t) 


(2.9) 


fi  2 

exp[-x  ✓ p“  + x + xp] 

2 2 
p + x 


f 

**  r\ 


U(t)  / JQ[x(t-u)J  JqIx  /ulu+2x)  ] du 
70 


(2.10) 


exp  (-T /plp+xj  xp]  <5(t) 


,x  [ 7 /*t(t+Tfy]  2 

♦ exp  (-=•  t)  — (~)  T 


j / t(t+2x)  ' 2 


(2.11) 


exp[-x  /p(p-fx)  + xp] 
p(p+x) 


U(t)  exp[-f]/  I0[f(t-u)]  I0[y/u(u+2x)]du  (2.12) 


J1(-ipx)  = - J^ipx)  = - i I^px) 


(2.13) 
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Jn(-ipx)  = Jn(ipx)  = In(px) 


Jn(-ipx)  = (-)n  Jn(ipx)  = exp(yiniT)  Wpx) 

i r271 

JQ(x)  =2 v J exp(ix  cos»)  d<J> 


J'(x)  = J , (x)  - - J (x)  = - J ,(x)  + -J  (x)  = 
n n-1  1 x n J n+lv  x nv 


7 - WX>1 


ft 

'a  w / 2Tt 
U P 


/X 

J^X)  dx  = -J- 
W T 


(The  substitution  uj  / 2rt  = x has  been  used.) 
P 


f JoK  /1™)du  " 4"  f 

JC\  * U>  T 0 


x JQ(x)dt  = 


x J:(x) 


, r xjocx) 

J W (u+2t)  )du  = — J ~7 

p pyo  /C 


X + 0)  T 


r r w (x  +<i>  T ) 

p 

(The  substitution  x = / t (t+2x)  and  (2.17)  have  been  used.) 

T hf  V“P /u('j*2T”du  - - T / * A ''^ldx  = 

yo  o L /77J? ] 

P 

t /**  <*>  T X 

■ - 1 Vx)  tTt  * v Vx)  ^ , , s/i  1 

o (x+uyr  ) 

(The  substitution  x = w / t(t+2x)  and  (2.17)  have  been  used.) 


Jn(x)  * /it  cos[x 


1 . TT  "I 

" * 2>  2 J 


for  x -►a®  and  n fixed. 


Legendre  and  associated  Legendre  polynomials 


(2.22) 


1'r’2  pLl"’>  ■ 2 " " '’Ll™  ’ 2 n p2ntn> 


P2„-l(0>  = 2 ” P2n“» 


(2.23) 


P2n(0)  = (-^  f0r 


(2.24) 


Spherical  Eankel  functions  of  the  second  kind 
(note  that  the  usual  superscript  (2)  has  been  omitted  for  notational  simplicity) 


h Cx)  = in  ix£(-ix)  ? 
nv  x L 


(n+m) ! 
m m! (n-m) ! 


(2ix)m 


(2.25) 


h2n-l(z)  = x h2n-2(z)  * h2n-2(z) 


(2.26) 


Modified  Lerah  function 


L(s,q)  = l n"^  exp(2?rins) 


(2.27) 


Upon  use  of  Mellin  transform 


L(s,q)  = 


L-q)expii(l-q)ir/ 


( 2tts) 


+ l -2— S--  „P('q~n-1  exp(iniT/ 2) 


(2.28) 


3 . Transient  radiation  from  elementary  electric  and  magnetic  sources 

Suppose  we  have  an  infinitesimal  oscillating  electric  dipole  in  a homoge- 
neous isotropic  environment,  whose  Fourier-transformed  moment  is  ue(oj) . The 
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dipole  is  located  at  the  origin  of  a spherical  coordinate  system  (r,0,<J>)  and  is 
oriented  parallel  to  the  direction  9=0.  The  components  of  the  dipole's  elec- 
tromagnetic transformed  field  are  given  by  the  well-known  expressions 


Er(u)  = C 


2 u (gj)  exp(-iojx) 


(3.1) 


G (»)  MpC-iut)  . n C(M2  iu  cl 

e9m  ■ 5 511,9  it  * r * -j) 

r 

ue(u)  exp(ia)T)  r(itt)  2 ia)1 

Vw)  3 4iF sin0  L“T—  + — J 


(3.2) 


(3.3) 


k k 

where  ? = (p/e)  , t =r/c  and  c = l/(ep)  . 

The  transient  fields  can  be  easily  obtained  for  the  case  of  a nondisper- 
sive  environment,  when  the  dipole  moment  varies  arbitrarily  in  time.  By  taking 
the  inverse  Fourier  transform  of  (3.1) 


- hf 


E (go)  exp(io)t)  du>. 


doing  the  same  for  the  other  field  components  and  letting 

/+  oo 

ue(u>)  exp(iwt)  dui, 

— OO 

it  is  easily  obtained 

rfi(f)  cU(t*)-i 
Er(t)  =^Fcose[-^ — J 

r [*  u (t*)  U (t*)  curt*)  1 

sin©l_  — s — + + — ^2 — J 


Mt) 


1 . JV**5  “e^l 

_sineL__+  __j 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


where  t*= t-r/c  = t-x  is  the  retarded  time,  and  a dot  means  differentiation  with 


respect  to  time. 
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The  corresponding  expressions  for  the  magnetic  dipole  are  immediately 


i 


k 


obtained  by  dulity  (H-*-E,  E -+-H,  £-*-l /C,  ug  -*■  - u^)  . Thus 


2 '“m1'*’! 

V'l  * 4i?C  cos9l.— * — T-J 

1 . JV'1)  cum(t*)  1 

H0(t)  * JiH  Sln9  [-T—  * ~ * — J 


v« 


4irr 


rurt*)  urt*)  I 

5in0[-V-  + -v-J 


(3.9) 


(3.10) 


(3.11) 


Assume  the  applied  signal  to  be  a pulse  of  time  duration  T.  Only  the  1/r 
terms  will  be  shown  to  be  important  at  distances  larger  than  cT  (see  Sec.  4). 
Furthermore,  a physically  realizable  dipole  can  be  assumed  elementary  as  long 
as  its  physical  dimensions  are  small  compared  to  cT. 

Let  C be  the  static  capacitance  and  L the  static  inductance  of  an  electric 
dipole  of  effective  height  h and  of  a loop  of  area  A,  respectively;  then  the 
relations 

(i)  dipole 

I(t)  = C V (t)  , ue  = - h I (3.12) 

(ii)  loop 

V(t)  = - L *I(t)  , um  = u A I (3.13) 

are  valid  in  the  limit  of  elementary  dipoles. 

Substituting  (3.12-3.13)  into  (3. 6-3.8)  and  (3.9-3.11)  respectively,  the 
result 

(i)  Electric  short  dipole  of  effective  height  h: 


EQ(t)  = C H (t) 


H (t)  = - sin0 

d>  J 4rrrc 


ChV (t*) 
4tttc 


sin0 


(3.14) 

(3. 15) 
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ttMF 


(ii)  Small  loop  of  area  A: 


hq  (t)  = 1 ^"'^2  sinQ 
4tttc 


m^sine 


(3.16) 


4tttc  L 


E*(t)  » - C Hq (t) 


(3.17) 


is  obtained  for  the  field  at  distances  larger  than  cT. 

Equations  (3.14-3.17)  clearly  show  that  the  radiation  fields  are  not  a 
faithful  replica  of  the  transmitted  signals  (currents  or  voltages  on  the 
radiators).  In  particular,  if  the  current  is  proportional  to  a Dirac  pulse 
6(t)  for  an  electric  dipole,  or  to  a step  function  U(t)  for  a magnetic  dipole, 
the  radiation  field  turns  out  to  be  proportional  to  S(t).  Therefore,  the  radia- 
tion field  flip-flops . 

4 . Radiated  power  and  radiated  fields  under  transient  conditions 

Under  steady-state  conditions  there  is  a clear  distinction  between  radiated 

and  reactive  fields.  Radiated  fields  are  the  1/r  components,  and  these  fields 

only  contribute  to  the  real  part  of  Poynting  vector  (in  frequency-domain)  in 

2 3 

lossless  media.  Reactive  fields  are  the  1/r  and  1/r  components  and  these 
fields  are  responsible  for  the  imaginary  part  of  the  Poynting  vector.  Accord- 
ingly, radiated  fields  are  associated  with  the  radiated  power  while  reactive 
fields  are  associated  with  the  imbalance  of  the  mean  stored  magnetic  and  elec- 
tric energies  within  the  period  of  the  oscillating  field. 

Let  us  now  examine  what  happens  under  transient  conditions. 

For  the  electric  dipole  the  Poynting  vector  £ has  r-  and  8-components 
given  by 

SrW’E9(t)H*(t)  5t[“2*F“u*T<F)2u2]}  '4a) 


r*  Z sin^9  cSfl’  c • 1 ,c  ,2  2l 

SeCt)  • - 6 (t)H  (t)  = - -J  - g JU  ♦ j-uu  . y(f)  u J 

(47TCr) 


(4.2) 
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where  the  subscript  "e"  has  been  dropped.  For  the  magnetic  dipole  the  expres- 
sions are  the  same  provided  that  the  magnetic  moment  is  considered  and  that  £ 
is  replaced  by  1/?. 

When  a time  integration  is  performed  between  t^  and  t2>  the  total  flux  of 
power  density  within  the  time  interval  (t^.t.,)  is  obtained.  This  flux  is  com- 
posed of  a term  proportional  to 

u2(t)  dt  (4.3) 

and  a second  term  depending  only  on  the  values  of  u,  u at  the  endpoints  of  the 
time-integration  interval.  This  second  term  can  always  be  made  equal  to  zero 
by  choosing  u(t^)  = uCt^);  u(t^)  = uCt^) . For  instance,  when  u(t)  is  identi- 
cally zero  outside  a finite  time  interval,  integration  between  (-®,  + °°)  makes 
this  term  vanish.  Accordingly,  this  part  is  identified  as  the  reversible  part 
of  the  power,  i.e.,  that  part  which  is  stored  in  the  medium  during  the  radia- 
tion process  but  can  be  recovered. 

On  the  contrary,  the  other  part  of  the  power  proportional  to  (4.3)  is  a 
non-decreasing  function  of  time.  Accordingly,  it  is  identified  as  the  irrevers- 
ible part  of  the  power,  or  radiated  part,  with  the  understanding  that  radiation 
means  an  irreversible  process:  the  radiated  power  is  absorbed  by  the  sphere  at 

infinity  and  cannot  be  recovered. 

Note  that  only  the  1/r  field  components  contribute  to  the  radiated  power. 
Accordingly,  they  will  be  labeled  radiated  fields.  The  other  fields  contribute 
to  the  reversible  part  of  the  power  and  will  be  labeled  local  fields. 

A further  characterization  of  the  fields  is  obtained  by  computing  in  which 
region  of  space  each  term  can  be  considered  as  the  dominant  one. 

Reference  is  made  to  a Gaussian  pulse  of  electric  moment 


P 

t. 


1 


T being  the  effective  width  of  the  pulse  (distance  between  inflection  points) . 


Substituting  (4.4)  in  (3.6)  one  obtains 


’ 2cu(t*)  f . cT  t*  /cT\21 

Er(t)  ■ Ts  “sS  L'“  ~ T-  * ( -)  J 

47IT  CT 


4nr  cT 


(4.5) 


E0W  - 4}  - 4 f f *(f)2]  <4-6> 


H*(t)  -f^2sin6[{(4f)2-4}-4f  f]  (4-7) 

4TTT  CT  ' 


It  follows  from  (4. 5-4.7)  that  the  radiation  terms  are  always  dominant 
provided  that 

r > > cT  (4.8) 

except  in  the  neighborhood  of  t*  = T/2,  where  the  radiation  terms  are  zero. 
Identical  results  are  obtained  for  the  case  of  a magnetic  dipole.  It  is  con- 
cluded that  the  spherical  surface  r = cT  divides  the  space  into  two  regions, 
which  can  be  labeled  near  space  (inner  region)  where  local  fields  are  dominant 
and  deep  space  (outer  region)  where  radiation  fields  are  dominant.  Accordingly, 
the  reversible  part  of  the  power  is  essentially  located  inside  the  near  space, 
while  the  irreversible  part  resides  and  propagates  in  the  deep  space. 


Let  us  now  consider  the  case  of  an  electric  dipole. 


Reference  is  made  to  Fig.  1.  Two  charges,  +q  and  -q,  are  separated  by  a 
distance  21.  At  a time  t = - Jl/v  the  charge  +q  is  suddenly  accelerated  and  is 
then  made  to  travel  at  constant  velocity  v,  so  that  it  collapses  on  the  charge 
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-q  at  time  t=£/v.  The  plot  of  the  electric  moment  u(t)  of  the  system  versus 
time  is  given  in  Fig.  1.  Accordingly,  a pulsed  current  is  flowing  in  the  time 


interval  |t  j < 1/ c,  from  z - Z to  z - - i. 


Fig.  1.  Relevant  to  physical  model  for  calculating  the 
transient  radiation  from  an  electric  dipole. 

Formal  expressions  of  current  density  in  the  time  and  frequency  domains 
are  the  following: 

J(r,t)  = - z qv5 (x)  S(y)  5(z+vt) 

* - z_  q<5(x)  <S(y)  <5(t  + ^)  , M < ^ 

(see  (2.1)) 

J(r,t)  - 0 , |t|>£  (5.1) 

J (r.oj)  = - z q<5  (x)  5 (y)  exp  F - iu:  7 J , | z | < Z 

J (r,oj)  * 0 , I z | > l (5.2) 


1 


z being  a unit  vector  in  the  z direction.  The  transformed  radiation  fields 
associated  with  (5.2)  are 


iAM  . it  , »W"t/<0An-cos6))  sin6 

9 v ' 2ttt  n n - cos© 


5 


n 


£ 

V 


(5.3) 

(5.4) 


By  inverse  Fourier  transforming  (5. 3-5. 4)  into  the  time-domain,  the  result 


EeCt)  ’ ' 4^  ^ole  ' [S(t*  * x)  ' 6Ct*  - X)1  (S'5) 

E 

H^(t)  = , x = 7 - §■  COS0  = | (n-cos0)  (5.6) 


is  obtained.  Equations  (5. 5-5. 6)  clearly  show  that  the  radiation  emanates  from 
the  end  points  of  the  configuration  at  times  -2,/v  (acceleration)  and  +2,/v 
(deceleration),  respectively  (see  the  diagram  of  Fig.  1).  Accordingly,  the 
radiated  field  is  not  proportional  to  the  time  derivative  of  the  current  (see 
(3.14-3.15))  when  a time  scale  less  than  22,/c  is  considered. 

The  transition  from  (5. 5-5. 6)  - microscopic  approach  - to  (3.14-3.15)  - 
macroscopic  approach  - is  readily  accomplished  by  letting  2£./c  approach  zero. 
One  obtains 


E0ct> 

d.(t) 


lim  c Sin0  - 

Q * 4tttc  2) 


21/ c -o 


5 E( 


g(t*-x)  . 

2x  ^ 


4^  sin06(t*)  (5.7) 


'9 


(5.8) 


thus  recovering  (3.14-3.15),  since  I (t) 2£  = - q6  (t)  2£  = - uQ5 (t)  . 

For  the  case  of  a magnetic  dipole,  reference  is  made  to  Fig.  2.  A con- 
stant current  Iq  is  suddenly  switched  on  at  time  t = - tg  and  switched  off  at 
time  t = +tQ.  The  current  flows  in  a circular  loop  of  radius  R located  in  the 
plane  (x,y). 
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Formal  expressions  for  the  current  in  the  time  and  frequency  domains  are 


the  following: 


I(t)  = In[U(ft.)  * U(t-tn)] 


siiuot. 


!(<*>)  * 2Iq  (5. 

The  transformed  radiation  fields  associated  with  (5.10)  are  those  of  a loop 
with  spatially  constant  current  (5.10).  These  are  known  to  be 


^ 0 

Hq(w)  = - sin(wt0)  J1  (^  R sin0)  exp(-iwT) 


%(“)  - * « H0 


where  J^(x)  is  the  Bessel  function  of  first  kind. 


observation  and  the  z axis.  No  radiation  emanates  from  points  BB ' (black  points) 
due  to  a cancellation  effect.  Accordingly,  the  radiated  field  is  not  propor- 
tional to  the  second  time  derivative  of  the  current  (see  3.16-3.17)  when  a time 
scale  less  than  R/c  is  considered. 


Fig.  2.  Relevant  to  physical  model  for  calculating  the 
transient  radiation  from  a magnetic  dipole. 


The  transition  from  (5.13-5.16)  - microscopic  approach  - to  (3.16-3.17)  - 
macroscopic  approach  - is  accomplished  by  letting  R/c  approach  zero. 

Reference  is  made  to  the  radiation  which  takes  place  at  t = - tg.  For  R-^O, 
it  seems  quite  reasonable  to  substitute  for  the  actual  radiated  waveform  some 
equivalent  pulses  (see  Fig.  3).  Reasonable  requirements  are: 

(i)  the  strength  of  each  pulse  must  be  equal  to  the  time  integral  of  the 
actual  waveform 
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(ii)  the  time  allocation  of  each  pulse  must  divide  the  actual  waveform 
into  parts  of  equal  area.  When  these  requirements  are  fulfilled. 


1 


(S. 14-5.17)  transform  into 


v t«ct*.v|x)-«(f.v!*>i 

HQ(t)»  lim  x — — 

9 R/c-0  2mC  X 


3 V* 


4tttc 


y sin06(t*+tQ) 


E,(t)  - - C HaCt) 


x * — sin0 
c 


(5.17) 

(5.18) 


thus  recovering  (3.16-3.17),  since  I(t)  =(J(t-t-tg),  apart  from  the  factor 
3/tt  * 0.995  = 1 (depending  upon  the  somewhat  arbitrary  identification  between 
actual  radiated  waveforms  and  equivalent  pulses) . 


Fig.  3.  Relevant  to  transition  from  continuous 
radiation  equivalent  pulses. 

6.  Heuristic  approach  to  transient  radiation  from  coaxial  apertures 

The  heuristic  approach  presented  under  Sec.  5 for  computing  transient 
radiation  will  now  be  applied  to  the  case  of  a coaxial  aperture. 

Reference  is  made  to  Fig.  4 where  the  radiation  from  an  open  ended  coaxial 
cable  of  circular  cross  section  is  considered.  However,  the  results  which  will 
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be  presented  can  readily  be  generalized  to  apply  to  radiating  apertures  of  dif- 


ferent geometry. 


Fig.  4.  Relevant  to  transient  radiation 
from  a coaxial  aperture. 


First,  the  macroscopic  approach  will  be  exploited.  Let  V(t)  be  the  inci- 
dent voltage  at  the  open  end  of  the  coaxial  cable.  The  electric  field  in  the 
aperture  has  only  a p-component  which,  under  the  assumption  of  perfect  open- 
circuit  conditions,  is  given  by 


Ep(t) 


ZV(t) 
p in  b/a 


C6.1) 


and  the  magnetic  field  is  zero. 

By  using  the  equivalence  theorem,  the  radiated  field  can  be  computed  as 
due  to  an  equivalent  magnetic  current  density  with  only  a ^-component, 

J .ft)  »E  (t) . A z-directed  electric  dipole  of  moment  given  bv 

m<j>  p * 


ue(t>  * IHtt / ? ™2 


V(t)  r^2  2, 

pe  InbTa  (b  ’a  5 


(6.2) 


2 2 

is  therefore  associated  with  the  aperture.  In  (6.2),  A = ir(b  - a ) is  the  aper- 
ture area  and  Z^  * [£  Jin  b/a] /2tt  is  the  characteristic  impedance  of  the  coaxial 
cable. 
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When  A is  small  compared  to  (cT)“,  the  higher-order  moments  associated 
with  the  aperture  are  negligible,  and  the  radiated  field  can  be  obtained  from 
the  results  of  Sec.  3.  Hence 


E0<«  ■ T, 


A_ 

2tt 


V(t*) 

4iTrc“ 


sin0 


E0(t) 

5 


(6.3) 


(6.4) 


Accordingly,  the  radiated  field  is  proportional  to  the  second  time- 
derivative  of  the  incident  voltage. 

As  a preliminary  to  the  study  of  the  radiation  from  the  microscopic  point 
of  view,  a model  similar  to  that  of  Sec.  3 for  the  electric  dipole  is  considered. 
A charge  q at  position  z = - l is  suddenly  accelerated  at  time  t = - z/v  and  is  then 
made  to  travel  at  constant  velocity  v along  the  positive  direction  of  the  z axis. 
At  t = 0 and  z = 0,  the  charge  is  suddenly  reflected  and  is  again  made  to  travel 
at  constant  velocity  v along  the  negative  direction  of  the  z axis.  At  t = z/v 
and  z = - l the  charge  is  suddenly  stopped. 

The  radiated  field  can  be  computed  by  following  the  same  procedure  as  in 
Sec.  5.  Only  the  results  are  given  in  the  following.  The  field  is  radiated  at 
times  t = - Jl/v,  0,  l/v,  from  points  z = -l,  0,  A,  respectively.  Particularly, 
the  field  radiated  at  t = 0,  z = 0 (reflection  of  the  charge)  is  given  by 


E0(t)  = ‘ 5 2?F 


nsin9 

2 2q 
n - cos  9 


6(t*) 


e9(t) 

c 


(6.5) 


(6.6) 


Turning  to  the  coaxial  structure  of  Fig.  4,  let  I (t)  =q6(t-z/c)  be  the 
current  propagating  along  the  cable.  The  pulsed  linear  charge  densities  on 
the  inner  and  outer  conductors  will  be  q/2TTa  and  -q2iTb,  respectively.  These 
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pulsed  charges  will  apparently  propagate  with  velocity  c along  the  cable. 
Accordingly,  the  radiated  field  can  be  computed  via  (6.5)  (with  n = 1)  and  by 
using  superposition.  Hence 

Ee(t)  ■ c - — -\  f 5(‘  - r>  - f s('  - r>l  (6-7> 

(2it)  r sine  IJQ  J 

where  r a ^ are  the  distances  of  the  point  (r,9)  from  the  inner  and  outer  rims 
of  the  cable,  respectively: 


r = r - a sin9  cosd> 
a 


r^  = r - b sin0  cos<(> 


(6.8) 


In  order  to  calculate  the  integrals  which  appear  in  (6.7),  it  is  convenient 
to  use  the  Fourier  transform 

/+<»  * 2tt 

dt  exp(-iojt)  / 5(t*  + — sin6  cos<t>)  d4> 

J0  c 


/+0°  b 

5(t*  + — sin6  cos<{>)  exp(-iwt)  dt 


= exp(-io)T) 


d<j>  exp(—  b sin6  cosifi) 


2 it  exp(-iu)t)J0(  ^ b sin6) , 


(6.8) 


and  similarly  for  the  other  integral.  Note  that  use  has  been  made  of  (2.16). 
By  Fourier  inverting  (6.8)  in  the  time-domain  and  substituting  in  (6.7),  the 


result 


Ee(t>  ' 4irtl*f  lEV>  - 1 


Mt) 


(6.9) 


(6.10) 


' 


is  obtained.  In  [6.9-6. 10)  the  functions  Ea’^(t*)  are  given  by 

Ex(t) 1—= j-r  , |t|  <|-  sine  (6.11) 

[(x  sinO/c)^  - t ]* 

Ex(t)  * 0 , |t|  > | sin6  (6.12) 

and  use  has  been  made  of  (2.4  and  2.14). 

A plot  of  the  radiated  field  is  given  in  Fig.  5.  The  preceding  results 
are  in  complete  agreement  with  those  of  an  exact  analysis  carried  out  by  using 
Wiener-Hopf  techniques  when  pulses  such  that  cT  < < b are  considered.  This  last 
restriction  does  not  apply  in  the  proximity  of  the  forward  direction. 


Fig.  5.  Plot  of  the  field  radiated  by  a coaxial  circular 
aperture  (sum  of  the  negative  and  positive  plots) 
and  transition  to  equivalent  pulses. 

Equations  (6.8-6.12)  clearly  show  that,  for  each  direction  of  observation, 
there  is  a continuous  radiation  for  |t*|  < b sin0/c  characterized  by  square- 
root  singularities.  The  radiation  emanates  mainly  from  points  AA',  BB'  ( flash 
points ) . Accordingly,  the  radiated  field  is  not  proportional  to  the  second 
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time-derivative  of  the  incident  voltage  V(t)  =Z^q6(t)  V6(t)  (see  (6. 3-6. 4)), 
when  a time  scale  less  than  b/c  is  considered. 

Also  in  this  case  the  transition  from  (6.8-6.12)  - microscopic  approach  - 
to  (6. 3-6. 4)  - macroscopic  approach  - can  be  accomplished  by  substituting  for 
actual  radiated  waveform  equivalent  pulses  (see  Fig.  5)  as  already  done  in 
Sec.  5.  Then  b/c  and  a/c  are  required  to  approach  zero.  One  obtains  (see 
Fig.  5) 


r V 

E„(t)  = lim  t i-r— 

9^  J . » , » Zn  4ttt  sin9 

b/c-*0;a/c-*0  0 


[(x-y)  ^ 


- x)  - 6 (t*  - y)  5(t*  + y)  - 3(t*  ♦ x) 

x-y  ' ^ yj  x-y 


. . C V b-a 

lim  — j—  

b/c-Oja/c-'-O  r /2c 


(x+y) 


3[t*  - ( (x+y) /2)  ] - <5[t*  + ( (x+y) /2) ; 


-i-  4-  ^sin9 

Z„  2-rr  4 2 

0 4irrc 


«*Ct)  ■ 


E0(t) 


b sinQ 


a sin9 


thus  recovering  (6. 3-6. 4),  since  V(t)  = V6(t). 


7.  Heuristic  approach  to  transient  radiated  from  wire  antennas 

Transient  radiation  from  wire  antennas  can  be  calculated  by  using  trans- 
mission line  techniques  for  modeling  currents  and  voltages. 

A cylindrical  antenna  of  height  2h  and  diameter  2a  is  considered  (see 
Fig.  6).  If  the  assumption  is  made  that  the  antenna  is  very  thin,  i.e.. 


= 2 £«[2h/a]  >>  1,  the  current  distribution  in  the  frequency  domain  reduces  to 


Fig.  6.  Plot  of  the  field  radiated  by  a wire  antenna 
fed  by  a rectangular  voltage  pulse. 


' = 2rr  q sin[oj(h-|z!)/c] 

A J in  -i  cos(>h/c) 


(7.1) 


where  Vin  is  the  transformed  voltage  at  the  antenna  terminals.  Equation  (7.1) 


is  formally  equivalent  to  that  of  the  current  distribution  along  an  open-ended 


transmission  line  of  length  h,  characteristic  impedance  ZQ  = 2tt  and  fed  by 


an  ideal  voltage  generator  of  zero  internal  impedance.  When  the  internal 
impedance  of  the  generator  is  different  from  zero  and  equal  to  aZ^  (real), 
(7.1)  can  be  immediately  generalized  by  using  the  preceding  transmission  line 
analogy.  Hence. 


Uz)  s 1 sin[h)(h-  1 z 1 ) /c] 

v Zq  a sin(ojh/c)-i  cos(ojh/c) 


(7.2) 


where  V is  the  transformed  voltage  of  the  generator. 

The  transformed  fields  associated  with  (7.2)  are  readily  seen  to  be 
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r - z 

y exp(-iu>r/c) 

cos  (wh/c) -cos  [ (ai/c)h  cos9] 

9 Z0(a+1) 

irr  sin0 

1+t1  exp(-2iuh/c) 

C 

(j  exp(-iwr/c) 

• (cos  (ioh/c) -cos  [ (w/c)  h cos9]} 

exp(-i  ^ h) 

' ZQ(l+a) 

rrr  sin9 

00 

• l (-nn 
on 

exp(-2niwh/c) 

(7.3) 

A 

F 

H - 9 

H ■ r 

(7.4) 

where  T = (l-ci)/(l+a)  is  the  reflection  coefficient  at  the  antenna  terminals. 
The  radiated  fields  in  the  time-domain  are  obtained  by  inverse  Fourier 

a 

transforming  (7. 3-7. 4).  For  a Dirac  pulse  excitation  V(t)  = V0(5(t),  V=Vq  and 
the  radiated  fields  are  given  by 


Ee^  - 


Z0(a+1)  27rr  sin9 


*{<5(t*)  + (1-r)  l (-Dn  6 (t*  - 2(n»1lh.) 

* „n  c 


- I C-r)n  [5(t*  - 


2n+l 


(h-h  cos9)) 


0 

+ 5(t* 


2n+l 


(h+h  co 


s6))]} 


E0(t) 


VtJ  - C 


(7.5) 

(7.6) 


Equations  (7. 5-7.6)  show  that  the  transient  radiation  of  the  antenna  con- 
sists of  a chain  of  Dirac  pulses,  each  being  a replica  of  the  applied  signal 
(save  for  the  amplitude  and  eventually  the  sign).  The  first  pulse  (first  term 
inside  the  curly  bracket  of  (7.5))  represents  the  radiation  due  to  the  launching 
of  the  current  along  the  antenna  at  time  t =0.  The  second  series  in  (7.5)  rep- 
resents the  radiation  from  the  upper  and  lower  tips  of  the  antenna  due  to  the 
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r 


charges  associated  with  the  current  and  going  back  and  forth  along  the  antenna. 

The  first  series  in  (7.5)  represents  the  radiation  due  to  the  currents 
relaunched  toward  the  generator.  These  currents  are  absorbed  in  the  internal 
resistance  of  the  generator  itself.  When  the  antenna  is  short-circuited  at  its 
terminals,  T=1  and  these  last  terms  disappear. 

For  an  arbitrary  applied  voltage  V (t) , the  radiated  field  is  easily 
obtained  from  (7. 5-7. 6)  by  using  convolution.  A sketch  of  the  transient  field 
radiated  by  the  antenna  when  the  applied  voltage  is  a rectangular  pulse  is 
given  in  Fig.  6. 

The  above  results  can  easily  be  applied  to  the  case  of  a monopole  antenna 
on  an  infinite  perfectly  conducting  ground  plane.  The  radiation  from  the  lower 
tip  of  the  antenna  disappears,  but  the  transient  pulses  reflected  from  the 
ground  plane  must  be  taken  into  account.  In  this  way,  a very  satisfactory 
agreement  with  experimental  results  is  obtained. 


8.  Heuristic  approach  for  transient  radiation  from  loop  antennas 

The  transient  radiation  from  a circular  loop  can  be  obtained  following 
similar  lines.  The  loop  of  radius  R is  fed  by  a voltage  generator  of  internal 
impedance  aZ^  (see  Fig.  7).  The  transformed  radiation  field  due  to  a traveling- 
wave  current  I^exp[-isx]  along  the  loop  is  known  to  be  [see  S.M.  Prasad  and 
B.N.  Das,  "A  circular  loop  antenna  with  travelling-wave  current  distribution," 
IEEE  Trans.  Antennas  Propagat.,  Vol.  AP-18,  No.  2,  pp.  228-280,  1970] 


VE1 


Jj(u>R  sin0/c)  oo  (i)  J^(u)R  sin9/c) 


2s 


L 


2 2 

s - n 


(s  cos  n<|)  - in  sin  rich) 


(8.1) 


oo  (i)  n J (ojR  sin0/c) 


VE1  COs6  £n  —2 27 

ln  (s  - n ) 


(ojR  sin0/c) 


(s  sin  n<j>  + in  cos  nd° 


r x E 


(8.2) 


(8.3) 
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[exp(-i2irs)  - 1]  exp[-iwx] 


Fig.  7.  Relevant  to  transient  radiation  from  circular  loop. 


Expressions  (8- 1-8.4)  are  the  generalization  of  the  field  radiated  by  a 
loop  with  a traveling-wave  current  Igexp[imx]»  m integer.  Letting  s=m+£  and 
then  5 approach  zero,  the  conventional  results  are  recovered. 

Now,  the  transformed  current  along  the  loop  is  assumed  to  be  coincident 
with  that  of  a transmission  line  of  length  irR,  characteristic  impedance  ZQ 
(real)  and  short-circuited  at  its  end.  Hence 

f 3 J cos[uR(TT-|xh/c] 

0 cos(wttR/c) 

V cos  [ajR(Tr-  lx!  )/c] 


q acos  (ojttR/c)  + i sin(ajiTR/c) 

A 

where  V is  the  transformed  voltage  of  the  generator. 

Expression  (8.5)  can  be  recast  in  terms  of  two  traveling-wave  currents 
moving  in  the  positive  and  negative  x angular  coordinate  directions.  By  alge 
braic  manipulation 


I = 


Va+1) 


exp(-iu)Rx/c)  + exp[-icaR(2iT  - x)/c] 
1 - r exp(-2Tiia)R/c) 


(8.6) 


From  (8.6)  and  (8. 1-8.4)  with  s-u>R/c,  and  using  superposition,  the  trans- 
formed field  radiated  by  the  loop  can  be  computed.  In  doing  that,  note  that 
for  the  second  traveling  wave  s =u>R/c  changes  into  -coR/c  and  E^  into  -E^.  Hence 


Jj(wR  sin9/c)  ® (i)n  J^(ujR  sin0/c) 


ojR/c 


♦ 2 - R • l 
c Jn 


(ooR/c)2  - n2 


cos  n<)> 


c _ p 2 cos9 
9 0 sin9  “n 


(i)  nJ  (wR  sin9/c) 

L 


(coR/c) 2 - n" 


sin  ncp 


r x E 


H = 


(8.7) 


(8.8) 


(8.9) 


where 
E0  = ' 1 


a 


2q  (ow- 1)  2toc 


[exp(-2iriujR/c)  -1]  expf-iojx]  £ Tn  exp(-2niriu)R/c)  (8.10) 


The  solution  for  the  transient  radiated  field  is  obtained  by  Fourier  transform- 
ing (8.10).  In  the  case  of  a Dirac  pulse  excitation  V(t)  =Vg5(t),  and  the 
transform  corresponds  to  the  field  radiated  by  a point  charge  going  around  in 
the  loop.  This  transform  can  be  easily  computed  by  using  results  (2.5,  2.7, 
2.8,  2.15,  2.17). 

9.  The  effect  of  dispersion 

Assume  now  that  the  electric  dipole  considered  in  Sec.  3 is  surrounded  by 
a homogeneous  isotropic  dispersive  environment  characterized  by  relative  per- 
mittivity £r((i>)  and  relative  permeability  pr(ai).  It  is  also  convenient  to 
introduce  the  relative  refraction  index  nr(u))  = /er((jo)yr(w)  . 

Letting  = /)Jq/£q  and  p = iw,  the  generalization  of  eqs.  (3. 1-3. 3)  to 
this  situation  is  the  following 
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2 uru(p) exp(-pnrt)  cos9 


r_E_  + — l 

rn  2 2 I * 
L r r n J 


U u(p)exp(-pn  t)  sin0  r 2 c *| 

Ee(p)  = S |_c  7ZT+  77  J 


yp) 
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(9.3) 


Assume  now  Ur  = l-  Relations  (9. 1-9.3)  can  be  recast  under  the  following 


| ex£(-£T)_sini9 

telpJ  4 nr 


[jrFlM 


where 


F : (p)  = u(p)  exp(-pt[nr(p)  - 1])  , 

exp(-px[n  (p)  - 1]) 


f2(p)  = U(P) 
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(9.4) 
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r 

(9.5) 
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(9.6) 

1])  , 

(9.7) 

1]) 

(9.8) 
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Eqs.  (9. 4-9. 6)  show  that  the  transient  radiation  from  the  dipole  can  be 
expressed  in  terms  of  the  transforms  of  the  two  functions  (9. 7-9. 8). 

Let  us  now  assume  u(p)  = Uq,  and  consider  the  case  of  a collisionless  cold 
plasma  for  which 


n (p)  = /l  - -f  • (9-9) 

P 

By  using  results  (2.9-2.10),  the  transient  radiated  field  can  be  cast 
under  the  following  form 
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’0  47TTC 


Uq  sin0 
’0  471TC 
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Uq  sin0 
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(9.12) 


where 


fx(t) 


2 J,(u)  /tTtZTr}) 

6 (t)  - uTt  — — p — U(t) 

p u)  /t(t+2f) 
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(9.13) 


f2(t) 


= U(t)  / J_[oj  (t-u)]  J [w  / u(u+2x) j du 


(9.14) 


It  is  noted  that  the  signal  is  identically  zero  for  t < t*  = t - T = t - r/c, 

which  is  the  usual  causality  condition. 

The  early  time  behavior  of  the  field  is  now  studied  at  large  distances 

from  the  source.  For  t-*-0,  w T -*-00,  the  function 

P 


J ^ ( co  / t(t+2x)  ) JjOjJ  / 2tt  ) 


u)  / t (t+2t) 
P 


w / 2tt 
P 


(9.15) 


tends  to  coalesce  toward  t = 0.  Use  of  (2.18)  shows  that  the  time  integral  from 

2 

-oo  up  to  +oo  equals  l/u^T.  Accordingly,  it  is  reasonable  to  substitute  to 

2 

(9.15)  a Dirac  pulse  of  amplitude  l/cy  and  centered  at  the  aenter  of  gravity 
of  the  area  of  (9.15)  (as  already  done  in  Sec.  5).  This  center  of  gravity  is 
readily  computed  from  (2.18)  to  be  equal  to  3/2.  Accordingly 


f (t)  * 6(t)  - 6 ft  - -V)  = -V  «(t) 

' 8oo  8u  T 


for  a)  T 
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Let  us  now  consider  the  function  f7(t).  For  t-*-0,  u t -*-00 


where  use  has  been  made  of  (2.19).  Again  the  function  J.(x)/x  approaches  a 


Dirac  pulse  as  u>  Accordingly 


Substituting  (9.16-9.17)  in  (9.10-9.12)  and  taking  into  account  (2.2),  the 


result 


Vin6  9 
0 4Trrc  „ 2 


is  obtained  for  t-*-0,  w t -*■<».  Accordingly,  only  the  radiative  components  of 
the  field  contribute  to  its  early  time  behavior,  and  the  field  tend  to  be 


vanishingly  small  for  oj 


The  late  time  behavior  of  the  field  is  now  studied  at  large  distances 


from  the  sources,  i.e.,  in  the  limit  u t*  -*■«>,  to  t -*-00,  t>>T.  From  (2.20-2.21) 


it  follows  that 


Substituting  in  (9.10-9.12)  and  using  (2.22)  the  result 
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is  obtained. 


(9.23) 
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Accordingly,  the  following  conclusions  are  drawn: 

(i)  The  magnetic  field  is  negligible  compared  to  the  electric  field. 

2 

(ii)  Both  the  1/r  and  1/r  components  play  a role  in  the  expression  of 
the  electric  field. 

(iii)  One  or  the  other  component  dominates  according  to  the  relative 

values  of  x/t  and  1/w  x. 

P 

The  case  of  a conductive  medium  with 


nt(p>  - /S 


(9.26) 


with  y -o/e^  and  a the  conductivity,  can  be  analyzed  in  a similar  way.  Tran- 
sient radiated  field  expressions  are  similar  to  (9.12)  with 


fL(t)  • «(t)  * (J)  x 


,2  Vi 


exp(J  t). 


(9.27) 


f2(t)  = exp£-  j t]  f (t-u)J  /u(u+2x)  ] . 


(9.28) 


10.  Transient  radiation  from  a spherical  antenna 

Let  us  consider  the  spherical  antenna  of  Fig.  8 referred  to  a spherical 
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system  of  coordinates  (r,  nscos9,  $)  centered  at  the  center  of  the  antenna. 

The  antenna  is  assumed  to  be  fed  by  a ^-independent  voltage  V (w)  (in  frequency- 
domain)  at  a small  equatorial  gap  of  thickness  2s,  s =a  aQ.  The  external  medium 
is  assumed  to  be  homogeneous,  isotropic,  non-dispersive  and  described  by  permit- 
tivity e and  permeability  u. 


I 


P(r,9) 


Fig.  8.  Geometry  of  the  spherical  antenna. 


The  boundary  value  problem  associated  with  the  radiation  from  the  antenna 
is  the  solution  of  steady-state,  source-free  Maxwell  equations  subject  to 

(i)  radiation  condition  at  infinity, 

(ii)  symmetry  condition  relative  to  the  (^-coordinate, 

(iii)  boundary  condition  at  r=a. 

The  solution  of  this  problem  can  be  obtained  following  the  same  procedure 
of  Ch.  1,3.  In  particular,  the  magnetic  field  has  only  a ({i-component  given  by 


4n-l 


~ ikVQ  oo 

Vr,9)  = ~T~  2 (2n-'l) 


h2n-l^ 


(2n-l)h2n_1(ka)  - ka  h2n_2(ka) 


P2n-1^  V2n-l<V 


(10.1) 


where  h^(x)  is  the  spherical  Hankel  function  of  the  second  kind,  k=u)v/"ql, 
S 3 / p/e  , P^(x)  the  associated  Legendre  polynomial  of  the  first  kind  and 
V2n  1^0^’  t*ie  v°ltaSe  excitation  for  mode  of  order  2n-l,  is  given  by 
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P2n-l(n)  dn: 


(10.2) 


this  is  a definition  slightly  different  from  that  of  Ch.  I,  eq.  (3.11).  Note 
that  for  a (spatially)  6-gap  excitation  vn=vn(°)  and  simplifies  upon  use  of 
(2.23)  as 


V2n-1  = P2n^ 


The  other  field  components,  i.e.,  Eg  and  E^.  immediately  follow  from 


(10.3) 


Maxwell's  equations: 


ikr^sinO  30  (sin9  V ’ 90  “ 


ife  h <*V 
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The  expression  of  the  field  in  time-domain  is  obtained,  in  principle,  by 
inverse  Fourier  (or  Laplace)  transform  of  (10.1  and  10.4)  once  the  Fourier  (or 
Laplace)  transform  of  the  applied  voltage  has  been  chosen. 

Assume  that  the  external  medium  is  free-space  and  that  the  applied  voltage 
changes  in  time  as  VgU(t),  i.e.,  as  a Heanside  step  function.  Then,  the  general 
term  of  the  series  expansion  (10.1)  that  must  be  inverted  is: 


h2n-l(“t/c> 


(2n-l)h2n_1  (wa/c)  - (aja/c)h2n_^(uja/c) 


(10.5) 


where  c = l//eQp0  . 

Firstly,  consider  the  term  n»l.  Use  of  (2.25)  transforms  (10.5)  into 
niia/c  + a/r 


a 

r 


iwa/c  + 1 - (o>a/c) 


/ . r-a\ 
2 exp(-iu)  — ) 


(10.6) 


which  exhibits  two  poles  in  the  complex  frequency  plane  p = ioi  at 


c . /3c 

21  1 1 SI  * - ai  ± lwl 


(10.7) 
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The  inverse-Laplace  transform  contour  of  integration  can  be  closed  in  the 

r.a 

left-hand  half  plane  for  t > and  the  integral  easily  computed  using  the 
residue  theorem.  Hence 


Hj>  (t)  = sine  expC-o^t*)  [cos(tt)jt*  «■£•)♦  sin  Wjt*] , 

t*  * t - ~ > 0 (10.8) 

At  large  distances  from  the  antenna  a/r-^0,  the  electric  field  has  only 
the  EQ1  = £ H component  and  the  radiated  field  is  a damped  oscillation  with 
period  4Tra/c  and  time  attenuation  constant  2a/c. 

Let  us  now  consider  the  general  term  (10.5).  Use  of  (2.25)  clearly  shows 
that  its  behavior  for  |p|  -*■<»  is  of  the  type  exp[p(r-a)/c]/p.  Accordingly,  the 
inverse-Laplace  contour  of  integration  can  be  closed  in  the  left-hand  half 
plane  for  t > (r-a)/c  and  the  inverse-Laplace  transform  computed  as  2ni  times 
the  sum  of  the  residues  at  the  poles  of  the  integrand.  These,  in  turn,  occur 
at  complex  frequencies  pg  given  by 

(2n-l)h2nl(-ipsa/c)  + i(psa/c)h2n_2(-ipsa/c)  = 0 (10.9) 

Note  that  the  (2n-l)  order  pole  of  h2n  ^(-ipr/c)  at  p = 0 cancels  with  the 
analogous  same  order  pole  of  the  denominator. 

Use  of  (2.25)  shows  eq.  (10.9)  to  have  2n  zeros;  each  residue  can  be 
easily  computed  upon  use  of  (2.26),  and  is  given  by: 


Ps  exp(pst)  h2n-l(*ipsr/c) 

2n(2n-l)  ♦ (p  a/c)2  h2n-l ("ipsa//c) 


U(t*), 


where 


t*  = t - 


is  the  retarded  time. 
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Relation  (10.10)  clearly  shows  that  the  dominant  poles  are  those  closest 
to  the  imaginary  axis.  On  the  other  hand,  the  analysis  of  Ch.  I,  6 showed  that 


these  poles  are  approximately  located  at 


P2n-Ia/C  = ' O^C211'1)273  1 i(2n-l) 


(10.12) 


for  n large.  Accordingly,  at  (retarded)  observation  times  large  compared  to 
(0.4) (2n-l) a/c,  the  first  few  modes  of  the  series  (10.1)  are  already  a 
good  representation  for  the  radiated  field.  Then,  the  conclusion  is  reached 
that  the  late-time  behavior  of  the  field  is  properly  described  by  the  inverse- 
Laplace  transform  of  the  first,  or  first  few  terras,  of  the  series  (10.1). 

Let  us  now  study  the  ear ly-time  behavior  of  the  field.  For  this  purpose, 
it  is  convenient  to  rearrange  the  general  term  (10.5)  to  be  inverted  under  the 
product  of  two  functions,  hence 


h2n-l(Wc) 


(2n-l)h2n_1(wa/c)  - (uia/c)  h2n_2(wa/c) 


h2n-lt“r/c) 
h.  " (ua/c) 


(10.13) 


The  inverse-Laplace  transform  of  the  first  term  can  be  immediately  obtained 
by  using  previous  results  as 


Ps  exp(pst) 
2n(2n-l)  * (psa/c)' 


(10.14) 


Analogously,  the  inverse-Laplace  transform  of  the  second  term  is  readily 
seen  to  be  equal  to 


- 6(t*)  - i - y 

r v J a Lc 


h2n-lC'ipqr/c) 

■h-nTn7-?/c)  u(t')’ 


(10.15) 


h2n  ^(-ip^a/c)  = 0,  upon  use  of  (2.25-2.26).  The  inverse-Laplace  transform  of 
(10.13)  is  now  obtained  by  convolution  of  (10.14  and  10.15).  However,  at  small 
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observation  times  only  the  convolution  with  the  6(t*)  function  will  be  impor- 
tant, so  that  the  field  radiated  by  the  antenna  will  be  given  by 


tr  • 


4n-l 

2(2n-l) 


P2n-1^ 


V2n-l^V  * Is 
1 


(Psa/c)  exp(pst*) 
2n(2n-l)  + (psa/c)2 


t*  -*-0 


+ 


(10.16) 


An  idea  of  the  convergence  properties  of  the  series  (10.16)  is  obtained  by 
considering  only  dominant  poles  (s  =2n-l),  the  field  in  the  equatorial  plane 
(H  = 0) , and  by  using  the  asymptotic  expression  of  the  series  term  for  n>>  1. 
Hence 


5Vn  r» 

H.  * =77  Re  £ 


exp[i(2n-l)ct*/a  - 0.4(2n-l) 


2/3 


ct*/a] 


(10.17) 


It  is  readily  seen  that  the  convergence  of  the  series  (10.17)  is  very  poor 
indeed  when  ct*/a  is  very  small,  and  that  the  series  diverges  for  t*  = 0. 

Considering  only  the  dominant  term  in  the  exponent  and  upon  use  of  (2.27), 
we  have  approximately 


5.8  V, 


H,  = 
<P 


TT£r 


1 


T7T 


(ct*/a) 


(10.18) 


Although  largely  approximated,  result  (10.18)  allows  us  to  conclude  that, 
at  small  observation  times,  the  inverse-Laplace  transform  of  the  series  expan- 
sion (10.1)  is  very  poorly  convergent  and  the  field  is  inversely  proportional 
to  a fractional  power  of  the  retarded  time  for  a step  time  excitation. 

On  the  other  hand,  the  early-time  behavior  of  the  field  can  be  computed  by 
using  results  of  Sec.  5.  As  a matter  of  fact,  at  time  t = 0 a magnetic  ring 
current 

Im(t)  = VQU(t)  (10.19) 
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is  applied  at  the  equatorial  plane  of  the  antenna.  For  t*  - 0 the  transient 
radiated  field  should  not  depend  on  the  particular  antenna  excited  by  the  mag- 
netic ring  current.  Accordingly  the  transient  field  should  be  coincident  with 
that  of  a magnetic  loop  on  a metal  surface.  By  modifying  results  of  Sec.  5 for 
9 =90°  and  by  using  duality  we  get,  for  ct*/a-*-0. 


/It r£r  /ct*/a 


(10.20) 


which  exhibits  behavior  of  the  type  (10.18).  This  proves  that  the  early  time 
radiation  from  the  antenna  essentially  comes  out  from  the  gap.  Using  the 
results  of  Sec.  S for  any  9 we  get 


0 1 


/2ttt  /ct*  sinS/a 


, ct*/a  sin0  < < 1 


(10.21) 


„ . _i  , t.  = t . 

4>  5 c 


(10.22) 


When  a/c  is  negligible  compared  to  the  time  resolution  of  the  measuring 
device,  the  procedure  developed  under  Sec.  5 transforms  (10.21-10.22)  in 


v 

En  = — sin9  6 (t*) 
0 TTrc  v ' 


(10.23) 


H . = -p— 

<P  S 


(10.24) 


since  only  one  half  of  the  magnetic  loop  contributes  to  the  radiated  field,  the 
second  half  being  screened  by  the  metal  body  of  the  antenna.  Accordingly,  the 
early-time  transient  field  is  proportional  to  the  time  derivative  of  the  applied 
voltage.  Incidently,  results  (10.21-10.22)  are  in  complete  agreement  with  the 
early-time  behavior  of  the  field  radiated  by  an  infinitely  long  cylindrical 
antenna. 


-63- 


CHAPTER  III 


; 


i 


TRANSFORM  PROPERTIES  OF  RADIATED  FIELDS 

Fourier , Fourier , how  many 
errors  are  made  in  your  name! 

I 

1 . The  problem 

In  elementary  antenna  theory  it  is  well  known  that  the  far  field  radiated 
from  a wire  antenna  is  just  the  Fourier  transform  - within  a proportionality 
factor  - of  the  current  distribution.  This  is  a rather  important  property  from 
the  applications  point  of  view.  It  is  helpful  from  the  analysis  point  of  view  - 
calculate  the  field  radiated  from  a given  antenna  current  - and  it  is  crucial 
from  the  synthesis  point  of  view  - from  a given  radiated  field,  calculate  the 
current  distribution. 

The  above  properties  of  radiated  fields  can  be  stated  in  an  elegant  way 
for  the  case  of  aperture  antennas.  It  is  assumed  that  an  aperture  exists  in  a 
metal  body  and  that  a known  field  distribution  is  produced  across  it.  When  the 
surface  of  the  metal  body  coincides  with  a coordinate  surface  of  a number  of 
coordinate  systems  - e.g.,  an  infinite  plane,  a circular  cylinder,  a circular 
cone,  a sphere  - it  is  possible  to  express  the  field  associated  with  the  aper- 
ture as  the  superposition  of  elementary  fields  constituting  a complete  basis. 

In  the  simplest  case  - an  aperture  in  an  infinite  metal  plane  - the  elementary 
fields  are  plane  waves,  so  that  the  plane-wave  expansion  of  the  radiated  field 
is  considered.  It  can  be  shown  that  this  plane-wave  expansion  allows  for  the 
directivity  diagram,  the  directivity  and  the  superdirectivity  ratio  to  be 
easily  computed. 

A further  transform  property  of  the  field  is  obtained  by  noting  that  the 
aperture  is  space  limited.  Accordingly,  it  can  be  shown  that  Hilbert- transform 
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relations  link  real  and  imaginary  parts  of  the  radiated  field.  This  is  a 
rather  important  property  from  the  practical  point  of  view,  since  it  reveals 
physical  constraints  to  which  the  radiated  field  must  comply.  Knowledge  of 
this  constraint  allows  for  realistic  synthesis  procedures  to  be  developed, 
for  measuring  procedures  to  be  simplified,  for  numerical  methods  and  computa- 
tions to  be  verified,  and  so  on. 

2.  Relevant  formulas  and  pertinent  expansions 
Transformations  of  "del"  operator 

7 • [A  exp(-i  Jt  • r)  ] * A • V exp(-i  k • r)  = -i  A • k exp(-i  k • r)7  • r = 

= -i  k_  • A expC-i  k_  • r) . Then,  formally 

7 • - -ik.  (2.1) 

7 x [A  exp(-i  • r)  ] = - Ax  V exp(~i  k_  • r)  = -ik  xA  exp(-i  Jc  • r) 

Then,  formally 

7 x -*■  -i  k x (2.2) 

Bessel  functions 

/oo 

exp(in£  + ixcos£)d£.  (2.3) 

. 

Asymptotic  evaluation  of  integrals 

Let  I(ft)  = f f(z)  exp[-iftq(z)]dz  with  q'(zs)  = 0,  q"(zs)  / 0;  zg  is  called  a 
G 

saddle  point.  Then,  in  the  neighborhood  of  z = zg, 

q"(z  ) 2 

q(z)  = qUs)  + — 2 — (2  - zs)  + — 

9Tbs)  2 2 

Letting  — * (z  - zg)  = s (stationary  phase  path) , for  ft  -*-  «> 
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provided  that  the  original  integration  contour  G can  be  deformed  into  the 


stationary  phase  path  and  the  saddle  point  is  sufficiently  far  from  singulari 


ties  of  the  integrand 


Similarly,  in  two-dimensions 


where 


Parseval 's  theorem  in  two- dimens ions 


A 

where  f(u,v)  and  g(u,v)  are  the  inverse  Fourier  transforms  of  f(x,y)  and 
g (x,y) , respectively. 


r . /to2  2 > 

1 fin  _____  \r  _ \r\r 

dv  - ^ rn-  S 

ay 

vb  > 

.2 

sin  XL_  c 

vb  \ 

sin(2b2  y > 

ay  - a cos  2 5 

L 2to 

< b 

TO  ' 

■“bill  -X  \-t 

2to" 

l b 

TO  ' 

where  C(x),  S(x)  are  the  Fresnel  integrals. 

Fresnel  integrals  and  auxiliary  funations 


2 2 

C(t)  = j + f(t)  sin(2|-)-  g(t)  cos  ) (2.10) 

2 2 

S (t)  =|  - f(t)  cos(H-j-  gCt)  sin(l|-)  (2.11) 

where  f(t)  and  g(t)  are  the  so-called  auxiliary  funations.  For  t-*0: 

f(t)  = j - J t2  + ...•  g(t)  = y - t + ...  (2.12) 

For  t 

f(t)  * ^ ; g(t)  = -3  (2.13) 

irt 

£’(t)  = - irt  g(t)  ; g’(t)  = irt  f(t)  - 1 (2.14) 


3.  Plane  wave  expansion  of  electromagnetic  waves 

Let  us  consider  the  electromagnetic  field  (E_,H) , solution  of  the  source- 
free  Maxwell  equations 


V x E = 

- i 60  p H 

(3.1) 

7 xH  = 

iueE 

(3.2) 

7 • E = 

0 

(3.3) 

7 • H = 

0 

(3.4) 

The  environment  is  assumed  homogeneous  and  isotropic,  and  steady-state 
conditions  are  postulated. 
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Assume  now  the  spatial  dependence  of  the  fields  to  be  of  the  type 
exp(-i  Jc*r).  Then,  upon  substitution  into  (3.1  through  3.4),  the  result 


k x E = co  y H 


(3.5) 


H x k = co  e E 


(3.6) 


k • E = 0 


k • H = 0 


(3.7) 


(3.8) 


is  obtained,  where  E and  H are  spatially  constant  vectors.  Note  that  use  has 
been  made  of  (2. 1-2. 2). 

Let  us  now  cross-multiply  eq.  (3.5)  by  k from  the  right,  and  substitute 
from  (3. 6) . We  get 


kxExk-w  £ y E = 0 


(3.9) 


~ 2 2 
Eq.  (3.7)  shows  that  E_  and  k_  are  orthogonal.  Then,  letting  w ey  = S , where 

6 is  the  free-space  propagation  constant. 


2 2 ~ 

(k  - 8 )E  = 0 


(3.10) 


Accordingly,  the  absolute  value  of  1c  is  determined,  and  equals  8-  Let  u,v,w 
be  the  Cartesian  components  of  1<: 


k_  = ux  + v£  + wz  ; 


(3.11) 


it  follows  from  k = 8 that 


772 72  27 

w = / 8 -(u  + v ) 


(3.12) 


Accordingly,  the  (u,v)  plane  can  be  divided  into  two  regions  separated 

by  a circle  centered  at  the  origin  and  of  radius  8<  For  the  inner  region, 

2 2 2 

which  is  called  the  visible  region,  u +v  <8  andw  is  real.  The  square-root 
in  (3.12)  is  taken  with  the  positive  sign,  so  that  propagation  along  the 
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positive  direction  of  the  z-axis  is  assured.  The  propagation  vector  k_  is 
determined  by  direction  cosines 


As  jc  components  are  not  all  independent,  also  E,H  components  are  mutually 


related.  In  particular 


(3.16) 


from  (3.5),  where  £ = /p/e  , and 


u E + v E 
x \ 


(3.17) 


from  (3.7). 

Under  rather  general  conditions,  an  electromagnetic  field  (E^H)  can  be 
expanded  in  a set  of  plane  waves.  Hence,  for  z>0 

J-H»  f +0O 

du  I dv  £(u,v)  exp(-iwz)  exp(-iux  - ivy)  (3.18) 

-UO  -00 

A 

Note  that  E_,  which  is  (u,v) -dependent,  is  a density  of  electric  field  per 
square-wavenumber , i.e.,  it  is  measured  in  volt*m.  Obviously,  the  formal  super- 
position (3.18)  is  of  practical  use  if  the  vector  function  E_  can  be  determined 
once  the  sources  of  IE  are  given. 

4.  Radiation  from  an  aperture.  Spectrum  calculation 

Let  us  now  cross-multiply  (3.18)  by  z_  from  both  left  and  right,  and  specify 
the  result  for  z = 0.  We  get 


J+00  M+  oo 

du  I dv  z x IB (u , v)  x z_  exp (-iux-ivy)  (4.1) 


-00  -00 


Note  that  z xE(x,y,0)  x.z_  = §t(x,y)  is  the  tangential  component  of  the  electric 
field  in  the  plane  z=0,  i.e.,  the  applied  tangential  electric  field  across 
the  aperture.  Letting  z xE(u,v)  x£  = §t(u,v),  eq.  (4.1)  can  be  rewritten  as 


► 4*00  M +00 


,v)  exp(-i  u x - ivy) 


(4.2) 


Eq.  (4.2)  shows  the  tangential  electric  field  across  the  aperture  to  be  the 
double  Fourier-transform  of  the  (x,y)  components  of  the  plane-wave  expansion 
of  the  electric  field  associated  with  the  aperture.  Accordingly,  by  inverting 
(4.2),  the  result 


Et(u,v) 


~+oo  -+00 

it?  j J dy 


y)  exp(i  u x + ivy) 


(4.3) 


is  obtained.  Note  that  the  integration  in  (4.3)  can  be  limited  to  the  aperture 
itself,  since  E£(x,y)  = 0 in  the  remaining  part  of  the  plane  z = 0.  The  vector 

A 

E£(u,v)  is  called  the  spectrum  or  spatial  spectral  components , associated  with 
the  aperture;  analogously,  k is  called  the  spatial  frequency  associated  with 
the  aperture.  The  third  component  of  the  complete  plane-wave  expansion  can 
easily  be  obtained  via  (3.17). 

If  E£(x,y)  can  be  factorized  into  two  functions  of  a single  variable,  i.e., 


E£(x,y)  = - Vt  fx(x)  f2(y). 


(4.4) 


(as  it  happens  in  rectangular  waveguides  for  TM  modes),  substitution  into  (4.3) 
and  integration  by  parts  leads  to  the  result 


§t(u,v)  = i Fx(u)  F2(v)kt 


(4.5) 


where 


i r00 

Fl(u)  = 2?  J exp(iux)dx 


(4.6) 


1 f +0° 

F2(v)  = 2tt  J f2  W exp(i  vy)dy 


(4.7) 


are  the  inverse  Fourier  transforms  of  f^ (x) , f2(y)  respectively,  and 
k£  = ux  + vy  is  the  transversal  component  of  the  wavenumber  vector. 
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Let  now  (p,40  be  polar  coordinates  in  the  z~0  plane,  and  let 
-t  = *^t  cos^*  + sin'^y-  Should  the  factorization 


Et(P,4>)  = - 7t  f(P)  exp (im}>) 


(4.8) 


I 


be  possible  (as  it  happens  in  circular  waveguides  for  TM  modes),  we  analogously 
get 


ik  r ^+27r 

§t(u,v)  = J ( f (p) P dp  | exp(in<)>)  exp(i  k.  • p) dc}>  = 

1 (2r 0 J J _t  - 


-ik  fCp)p 

(2ttj 


f f(p)pdp  f 

0 0 


2ir 


dp  I exp(inC  + i ktp  cos^)d^  (4.9) 


Then,  by  using  (2.3),  the  result 


§t  (u,v)  = iktir 


exp  (ini);) 


if 


f(p) 


Jn(ktp)p  dp 


(4.10) 


is  obtained.  Accordingly,  E is  factorized  in  terms  of  two  functions  of  ip 
and  p respectively,  the  latter  being  the  inverse  Fourier-Bessel  transform  of 
the  aperture  field  function  distribution  f(p). 


5.  Radiation  from  an  aperture  continued.  Far-field  calculation 

Let  us  investigate  how  the  integral  (3.18)  simplifies  when  the  electric 
field  is  evaluated  at  large  distances  _r  from  the  aperture.  A particular 
spatial  direction  is  determined  by  its  directional  cosines  (£, n,a): 

r = r(£x  + n£  + oz)  (5.1) 

2 2 2 

where  £ + n + o =1.  Along  this  particular  direction,  at  a large  distance  from 
the  aperture,  the  electric  field  Il(r)  may  be  written  as  the  superposition  of 
elementary  plane  wave  fields  according  to  (3.18).  However,  at  large  distances 
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from  the  aperture,  the  field  is  a locally  plane  wave.  Accordingly,  we  can 
reasonably  surmise  that  the  field  will  be  mainly  set  up  by  those  plane  wave 
constituents  whose  propagation  vectors  are  contained  within  a small  cone  whose 
axis  is  the  considered  direction. 

In  order  to  make  this  intuitive  reasoning  mathematically  sound,  let  us 
compute  the  value  of  the  argument  of  the  exponential  in  (3.18),  namely 
g(u,v)  = ux  + vy  + wz,  when  the  propagation  vector  1c  is  oriented  in  the  neigh- 
borhood of  the  considered  direction  (C,n,o).  This,  in  turn,  requires  the  first 
two  cosine  directions  of  k to  be  in  the  neighborhood  of  u^.v^,  where 

uQ  = SC  5 vQ  = 6n  (5.2) 


Accordingly,  we  have 

g(Vv0)  = Br  (C2  + n2+a2)  = Sr 


(5.3) 


aC 


A 


n / 


-v° 


(5.4) 


(5.5) 


It  follows  that  g(u,v)  is  stationary  in  the  neighborhood  of  (u  ,Vq).  The  main 
contribution  to  the  integral  (3.18)  then  arises  from  an  integration  area 
centered  at  (Uq,Vq)  and  vanishingly  small  as  Sr-*-°°.  As  a matter  of  fact,  out- 
side this  area  the  phase  of  the  integrand  is  a rapidly  varying  function  of 
(u,v)  . Accordingly,  as  Sr-*-00 


E(r,9,cf>) 


— (Uq  « vQ) 


exp(-i  Sr) 


dv  exp [-i (1c  • £ - Sr) ] = 


S exp(-i  Sr) 


cos9  E(u0,vQ) 


2ni 


r 


(5.6) 


Note  that  use  has  been  made  of  (2.5)  for  evaluating  the  integral 


Relation  (5.6)  states  that  the  far- field  radiation  diagram  of  the  aperture 


is  determined  - aside  from  angular  factor  cos9  - by  its  plane-wave  associated 


expansion.  For  a given  (spherical  coordinates)  direction  (9 , <+>)  the  values 
(Uj-.Vp.)  appearing  in  (5.6)  are  given  by 


Let  us  now  consider  a rectangular  aperture  of  dimensions  (2a, 2b).  The 


is  introduced,  centered  at  the  center  of  the  aperture,  with  the  x-axis  parallel 


to  the  side  of  length  2a  and  the  y-axis  parallel  to  the  side  2b.  The  aperture 


is  illuminated  by  an  electromagnetic  field  whose  tangential  electric  field  com 


ponent  is  given  by 


The  tangential  component  E of  the  spectrum  assoicated  with  the  aperture 


is  readily  computed  as 


the  z-component  of  the  spectrum  follows  from  (3.17),  hence 


Let  us  now  study  the  radiation  diagram  in  two  principal  planes 


The  E-plane  is  defined  by  d>  = - 90°.  For  this  plane  we  have,  from  (5.7), 


uQ  = 0,  vQ  * ± g sin0.  Then,  from  (6.3) 


E = ± tanQ  E 

z y 


Accordingly,  the  far- field  will  have  components 


E0  = ' i E0 


H$  = 5 


8 ab  , .0  , fsin(0  b sin0)"l 

^ exp('l6r)  L s.sue  J 


(6.4) 


(6.5) 

(6.6) 


where  use  has  been  made  of  (5.6). 

The  bracketed  term  in  (6.5)  is  clearly  the  radiation  diagram  in  the  E-plane. 
The  first  radiation  null  occurs  at  an  angle  0^  given  by 

6 b sin0Q  = 7 r.  (6.7) 

For  large  apertures  (0b>>l),  0^  is  small  and  the  width  of  the  first  lobe 
is  given  by 

% - F (6-8) 

The  first  lateral  lobe  occurs  at 

0b  sin0  = (6.9) 

Its  amplitude  is  equal  to  2/377,  i.e.,  lateral  lobes  are  more  than  14db  below 
the  principal  lobe. 

The  H-plane  is  defined  by  <j>  = 0,  180°.  For  this  plane  we  have,  from  (5.7), 

Uq  = ± B sin0,  Vq = 0.  Then,  from  (6.3),  E2  = 0 and  the  far-field  will  have  com- 
ponents 
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i Eq  |-^r  cos6  exp(-iBr)  — 


cos  (B  a sin8)  ~ 
/ 2 B a sin8  \- 


(6.10) 


where  use  has  been  made  of  (5.6). 


(6.11) 


The  bracketed  term  in  (6.10)  times  cos8  is  clearly  the  radiation  diagram 
in  the  H-plane.  The  first  radiation  null  occurs  at  an  angle  80  given  by 


B a sin8g  = y- 


(6.12) 


Note  that,  at  B a sin8  = tt/2,  the  radiation  diagram  is  not  zero,  taking  instead 
the  value  -rr/4 . For  large  apertures  we  have 


AO  ~ i! 
A0O  - 2a 


(6.13) 


The  first  lateral  lobe  occurs  at 


B a sin8  = 2ir 


(6.14) 


Its  amplitude  is  equal  to  1/15,  i.e.,  lateral  lobes  are  more  than  25  db  below 
the  principal  lobe. 

Previous  discussion  highlights  the  fact  that  an  amplitude  tapering  of  the 
field  in  the  aperture  reduces  the  amplitude  of  lateral  lobes  at  the  expense  of 
an  increase  in  the  width  of  the  principal  lobe. 


7.  Example  continued.  The  effect  of  phase  tapering 


The  effect  of  phase  tapering  in  the  aperture  is  now  examined.  To  this  end, 
the  tangential  component  of  the  electric  field  in  the  aperture  is  assumed 


equal  to 


Ea  = E, 


y sin(i  x)exp[i  ^-(i  - ri) 
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Comparison  of  (7.1)  with  (6.1)  shows  that  the  u-dependence  of  Eg(u,v)  is 
the  same  as  in  (6.2),  while  the  v-dependence  is  given  by 

b 2 

Q(v,a)  = exp(i  j a2)  J exp^i  vy  - i y a2  ^ ) dy  = 


= l a)-  1 f(^-  °)]  exP(ivb) 

-[4S+  a)‘  1 f(S+  a)]  exp(-i  v b)| 


(7.2) 


where  the  functions  f and  g are  defined  by  (2.10-2  11).  If  the  further  assump- 
tion that  a is  small  is  made,  the  functions  f and  g can  be  expanded  in  the 
neighborhood  of  \\b/ira  and  only  the  first  two  terms  of  the  expansion  retained. 
Then,  upon  use  of  (2.14),  the  result 


«*.«>  *(S)*  !?*(£)] 


(7.3) 


is  obtained. 


The  effect  of  the  (quadratic)  phase  tapering  in  the  E -plane  is  now  studied. 
Two  angle  ranges  are  considered.  Firstly,  assume 


8 b sin6  _ 2b  sin9 


(7.4) 


i.e.,  angles  0 much  smaller  than  a.  Then,  upon  use  of  (2.12), 


Q(v,a)  - 2b  cos(8bsin9) 


(7.5) 


Accordingly,  the  radiation  diagram  is  not  appreciably  changed  with  respect  to 
the  case  of  uniform  phase  illumination. 

Conversely,  assume 


2b  sin9 


» 1, 


(7.6) 
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then,  using  (2.133,  expression  (7.3)  reduces  to 


Q(v,a)  = 2b 


fsin  v b 

. / TO  1 

L v b 

H^b)  cosvbJ 

(7.7) 


Accordingly,  the  radiation  nuls  are  "filled  up". 

The  previous  discussion  shows  that  a quadratic  phase  tapering  (or  quadratic 
phase  error,  as  is  usually  called)  tends  to  smooth  the  radiation  diagram. 


8.  Directivity  and  superdirectivity  ratio 

The  (complex)  power  associated  with  the  aperture  is  given  by 


i c c i r°°  f+°° 

P = 2 dx  I dy  ExH*  • | = y dx  I dy(z  xE)  • H * = 

-00  -CO  •'_£» 


loo  •'-oo 


*r*r 


dy(z  x E)  • (z  xH*  x 2) 


(8.1) 


oo  -00 


By  applying  Parsevals’  theorem  and  (3.5),  eq.  (8.1)  can  be  transformed  as 
follows 


+00  r +00 


•u 


dy(zxE)  x (zxH)*  • z = 


00  -00 


* 2^1  J du  I dv(2x§)  x (£xkx|)**z 
— 00  —00 


(8.2) 


Relation  (8.2)  can  be  recast  in  a different  way  by  letting  E = E + E z, 

~ 1 ~ t 2 ~ 

A 

k = + w z,  evaluating  the  vector  products  and  recalling  that  u and  v are 

real.  Then 


2 f +00  .+00  _ 

P=1Sr/  du|  dv[w*  |lt|2  + w(Ej2] 


(8.3) 


-00  <-00 


Let  us  now  distinguish,  in  (8.3),  the  integration  in  those  regions  for 
which  w is  real  ( visible  domain ) and  imaginary  ( invisible  domain),  respectively. 


Calling  CR  and  Cj 

P - mil 

- 2? 


the  two  domains  respectively,  eq.  (8.3)  transforms  into 


(8.4) 


The  first  integral  which  appears  in  (8.4)  is  clearly  the  (real)  radiated 
(or  active)  power  from  the  aperture;  the  second  integral  is  the  reactive  power 
associated  with  the  aperture,  i.e. , l/2w  times  the  difference  between  the  time- 
averaged  values  of  the  magnetic  and  electric  type  energies  stored  in  the  half- 
space z > 0.  Accordingly,  that  part  of  the  integral  (6.4)  containing  |§t|^ 

■ a ,2 

represents  the  (magnetic  type)  stored  power  while  that  part  containing  |E  | 
represents  the  (electric  type)  stored  power.  The  total  stored  power  is  given  by 


(8.5) 


The  directivity  ratio  of  the  aperture  is  defined  as  the  (spatial)  density 
power  radiated  in  the  direction  of  main  radiation  divided  by  the  (spatial)  mean 
radiated  power.  Assume  the  main  radiation  to  take  place  along  9 =0°,  i.e.,  in 
the  direction  of  the  z-axis.  Then,  using  (5.6  and  8.4),  we  get 


4tt  g3  | £ (0, 0)  |2 
ff  w |§(u,v) | ^ du  dv 


(8.6) 


Note  that  the  (spatial)  mean  radiated  power  has  been  computed  with  reference  to 
the  solid  angle  4u. 

The  superdirectivity  ratio  y of  the  aperture  is  defined  as  the  total 
(radiated  plus  stored)  power  associated  with  the  aperture,  divided  by  the 
radiated  power.  Using  (8. 4-8. 5)  we  get 
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Assume  now  the  value  ji(0,0)  to  be  fixed  so  that  the  value  of  the  radiated 
field  in  the  direction  0 = 0 is  determined  via  (5.6).  Note  that  the  electric 
field  is  normal  to  z along  this  particular  direction.  From  (4.3): 

|(0,0)  = E (0,0)  = — ^ f dx  f dy  E*(x,y)  (8.8) 

r (2ir)  J J z 

so  that  the  (spatial)  mean  value  of  the  applied  electric  field  across  the  aper- 
ture is  fixed. 

In  order  to  increase  the  directivity  D,  the  integral  which  appears  in  the 
denominator  of  (8.6)  should  be  kept  as  low  as  possible.  This  implies  that  the 

A 

spectral  components  of  E_  in  the  domain  should  be  kept  as  low  as  possible. 
This  result  can  be  accomplished  by  illuminating  the  aperture  with  a tangential 
field  E_  which  exhibits  a (spatial)  rapid  variation,  so  that  the  (spectral) 
low-frequency  components  of  its  associated  spectrum  will  be  lowered.  Note  that 
the  (spatial)  mean  value  of  E3  should  be  different  from  zero. 

This  type  of  aperture  illumination  has  the  consequence  that  the  (spectral) 
high-frequency  components  of  the  associated  spectrum  will  be  enhanced.  Accord- 
ingly, the  superdirectivity  ratio  will  be  increased  and,  consequently,  the 
*3red  electric  and  magnetic  power  in  the  neighborhood  of  the  aperture  will 
• <*n  very  sensitive  to  a change  in  the  (temporal)  frequency  of  the  applied 
Vs  a consequence,  the  transmitted  bandwidth  will  be  decreased.  In 
the  superdirectivity  ratio  should  be  chosen  according  to  the  band- 
- • in  emitted;  this  results  in  an  upper  limit  to  the  possible  obtain- 
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9.  Hilbert  transform  properties  of  the  fields 


Let  f(x)  be  a square-integrable  function  such  that  f(x)  =0  for  x<0  and  let 


f(u) 


_1_ 

2tt 


exp( i u x)dx 


(9.1) 


A 

be  its  inverse  Fourier  transform.  Let  z=u  + iu'and  f(z)  be  the  analytic  con- 

A 

tinuation  of  f(u)  onto  the  conplex  E-plane.  Inverting  (9.1)  we  get 


f(x) 


= J £(u) 


exp(-i  ux)dx 


(9.2) 


The  integration  contour  which  appears  in  (9.2)  can  be  conveniently  closed  by 
means  of  a circle  of  infinite  radius  in  the  half-plane  u'  > 0,  for  x<0.  Then, 
since  f(x)  =0  for  x < 0,  it  follows  that  f(z)  should  be  singularity-free  in  the 
half-plane  u'  > 0,  (the  real  axis  is  included). 

Let  us  now  consider  the  integral 


J 


f(z) 

z-u„ 


dz 


(9.3) 


along  the  contour  C depicted  in  fig.  2.  Since,  from  (9.1),  f(z)  tends  uniformly 
to  zero  as  |z|  -+00,  the  integration  contour  C can  be  conveniently  closed  by 
means  of  a circle  of  infinite  radius  in  the  half-plane  u'  >0.  However,  f(z) 
is  singularity-free  in  this  half-plane.  Accordingly,  the  integration  along  the 
above  closed  contour  equals  zero,  i.e., 


ii hi 

u-u„ 


- TT  i f(U-) 


= 0 


(9.4) 


where  the  integral  is  to  be  understood  as  the  principal  value  and  the  second 
term  in  (9.4)  accounts  for  one-half  of  the  pole  contribution  at  u = Uq.  Now, 

A /v  -A 

if  real  and  imaginary  parts  of  f(u)  are  separated,  namely  if  f(u)  = f (u)  +i  f.,(u), 
relation  (9.4)  splits  into 


Fig.  2 - Relevant  to  the  Hilbert  transform  definition 

The  generalization  of  (9.4)  to  a two-variable  function  f(u,v),  double 
inverse  Fourier  transform  of  f(x,y)  such  that  f(x,y)  =0  for  any  x<0  and  any 


Relations  (9. 5-9. 6)  link  real  and  imaginary  parts  of  a function  which  is 
analytical  in  the  half-plane  u'  < 0 and  tends  uniformly  to  zero  as  |u|  -*-00.  How- 
ever, their  use  is  rather  impractical  if,  let  us  say,  the  imaginary  part  should 


be  recovered  from  the  knowledge  of  the  real  part.  Accordingly,  it  would  be 
helpful  to  obtain  approximate  relations  between  real  and  imaginary  parts  at 
the  same  abscissa  un.  To  this  end,  let  us  consider  an  integral  of  the  type 


/•  +00  r 0 

i = £ sM.  du  a f zM 

J u-u«  J u-u. 


(9.10) 


where  the  function  g(u)  is  explicitly  assumed  to  be  even  and,  without  loss  of 
generality,  Uq  has  been  taken  positive. 

Now  let 


u = un  exp(t)  , u > 0 


(9.11) 


u = - un  exp(t) , u < 0 


(9.12) 


Substituting  (9.11-9.12)  into  (9.10),  splitting  the  integration  path  of  the 
second  integral  in  (-°°,Uq)  , (Uq,oo)>  integrating  by  part  and  rearranging  terms 


results  in 


,+0°  dg(-uQ  exp[t]) 


in  [l  + exp(t) ]dt 


■+0°  dg[u  exp  [t] ) 


in 1 1 - exp(t) [ dt 


in  coth  — dt 


(9.13) 


Result  (9.13)  is  still  exact  under  the  only  hypothesis  that  g(u)  is  even.  Note 
also  that  the  sign  of  principal  value  can  be  omitted,  since  the  integral  of  the 


function  In  coth  |t|/2  from  to  +“>  is  finite  and  equal  to  tt2/ 2 . This  func- 
tion peaks  rather  sharply  in  the  neighborhood  of  t = 0 and  can  be  approximated 
by  a Dirac  pulse  of  area  tt2/2.  Then, 


It  is  noted  that  E^  is  just  the  inverse  Fourier  transform  of  the  tangent  elec- 
tric field  in  the  aperture.  Accordingly,  §t  is  singularity-free  in  the  half 
space  u'  >0,  v'  >0.  The  two  factors  k^  and  w do  not  introduce  any  further 
singularity  so  that  expression  (9.16)  still  shares  this  property.  However,  the 

A 

asymptotic  behavior  of  (9.16)  is  different  compared  to  that  of  E^,  so  that,  in 
general,  real  and  imaginary  parts  of  the  radiation  diagram  are  not  a pair  of 
Hilbert  transforms.  However,  §t(u,v)  decays  more  rapidly  than  l//u~  + v^  for 
| u | and/or  |v|  -*00  if  a spatial  tapering  in  the  aperture  is  present.  In 
this  case,  reasonings  leading  to  Hilbert-transform  properties  of  analytical  func- 
tions apply,  and  real  and  imaginary  parts  of  the  radiation  diagram  are  a Hilbert- 
transform  pair. 
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APPENDIX 


I . 

Exercises 


1.  Derive  Hallen's  integral  equation  for  a thin  wire  antenna  fed  at  its  cen- 
tral point  by  an  ideal  current  generator. 

Hint:  Express  the  ^-component  of  magnetic  field  in  terms  of  the  vector 

potential.  Substitute  the  expression  of  the  vector  potential  in  terms 
of  the  current  distribution  along  the  antenna.  Then,  specialize  the 
result  for  z = 0 and  take  the  source  specification  into  account. 

2.  Consider  an  infinitely  thin  current  filament  extending  along  the  z-axis 
from  -<»  to  +».  The  amplitude  I of  the  current  is  constant  along  z and 
varies  in  time  as  exp(jo)t).  Derive  the  electromagnetic  fields  associated 
with  this  idealized  antenna. 

Hint : Solve  Maxwell’s  equations  for  the  magnetic  field  and  remember  that 

(x)  = as  x+0. 

1 J -rrx 


3. 


Discuss  the  transient  current  response  of  a wire  antenna  of  length  2a  and 
diameter  2b  driven  by  a gaussian  voltage  pulse  v(t)  = Vgexp [- 2t^/T“] . The 


pulse  "width"  T is  larger  than  m/£Q|jQa.  The  antenna  is  surrounded  by  a 
homogeneous  isotropic  collisionless  cold  plasma  whose  relative  permittivity 

is  given  by  £(;<//)  =1  - . 

ur 


The  antenna  is  surrounded  by  a 


is  given  by  £(jw)  =1 


Detailed  computations  are  not  required,  although  a clear  statement  of  the 
problem  is  expected. 


Hint:  Use  the  equivalent  circuit  representation  for  the  input  admittance 

of  the  antenna,  truncated  after  the  first  dynamical  term.  For  a better 


understanding  of  the  branch-cut  singularity,  assume  a small  collision  rate 
v in  the  plasma,  so  that  £(jw)  =1  - . 


4.  Consider  an  infinitely  long  cylindrical  antenna  of  radius  p=a.  The 

antenna  is  assumed  perfectly  conducting  and  driven  at  z = 0 by  a voltage 
V exp(jiot)  applied  to  a 6-gap.  Derive  the  modal  expansion  for  the  mag- 
netic field  generated  by  the  antenna.  Give  the  formal  expression  of  the 
input  admittance  of  the  antenna. 

Hint : Assume  the  <j>-component  of  magnetic  field  to  have  a z-dependence 

exp(-ju|z|),  0<u<®.  Using  Maxwell's  equations,  verify  that  the  radial 

dependence  is  of  type  h|^ ( /k^  - u^  p) . Use  superposition  to  obtain  the 
general  expression  for  the  magnetic  field. 
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5.  An  elementary  electric  dipole  of  length  h and  capacitance  C is  driven  by 
a current  IgU(t),  where  U(t)  is  the  Heaviside  step  function.  At  large 
distance  r there  is  a receiving  elementary  dipole  of  length  h and  capaci- 
tance C,  parallel  to  the  transmitting  one.  The  receiving  dipole  is 
closed  on  a purely  conducting  load  G.  Compute  the  voltage  across  the 
load  as  a function  of  time. 

Hint:  Neglect  radiation  conductance  of  the  dipoles.  Take  into  account 

only  the  radiation  term  in  the  transmitted  field. 


6.  The  elementary  electric  dipole  of  figure  is  driven  by  triangular  current 
pulses  such  that  cT>>h.  The  dipole  is  vertically  located  and  at  a dis- 
tance d from  a perfectly  conducting  plane.  Find  the  direction  of 
maximum  and  minimum  (zero)  radiation. 


♦ I(t) 


Hint : Apply  the  image  theorem. 


7.  A wire  antenna  of  length  2a  and  diameter  2b,  a»b,  is  driven  by  a voltage 
pulse  v(t)  = Vq  f(t) , f(t)  = 1 forO<t<T,  f(t)  = 0 otherwise.  Assume  cT 
large  compared  with  b.  Calculate  the  input  current  as  a function  of  time. 

Hint:  Use  the  equivalent  circuit  truncated  after  the  first  dynamical  term. 


8.  The  elementary  electric  dipole  of  figure  is  driven  by  a current  I (t) . A 
passive  reflector  is  located  at  a distance  d from  the  dipole.  Assume  the 
reflector  to  be  in  the  near  field  of  the  transmitting  dipole.  Compute 
the  field  radiated  along  the  end-fire  direction  z > 0. 


Hint:  Consider  the  reflector  as  a receiving  antenna  closed  on  a short- 

circuit. 


r 


9.  A loop  antenna,  of  radius  R and  inductance  L,  is  located  parallel  to  a 
perfectly  conducting  plane  and  at  a height  d above  it.  The  loop  is 
driven  by  a guassian  pulsed  voltage 

V(t)  = VQ  exp(-2t2/T") , such  that  (cT)2»7TR2. 

Calculate  the  transient  far  field. 

Hint:  Apply  the  image  theorem. 

10.  Consider  a rectangular  aperture  in  an  infinite  conducting  plane  (see  figure) 


The  aperture  is  illuminated  by  an  electromagnetic  field  (£  ,Ha)  of  angular 
frequency  oj,  and  the  tangential  electric  field  distribution  across  the 
aperture  is  given  by 


Ea  = Eq  cos(^J)  exp(ihgy)  y. 


Determine  the  value  of  h for  which  the  radiation  diagram  has  a maximum  at 
0=30°  in  the  E-plane  (0  is  the  angle  with  respect  to  the  z-axis,  and  the 
E-plane  is  he  x = 0 plane) . 

Consider  the  rectangular  aperture  of  exercise  #10.  Determine  the  dimen- 
sion 2b  such  that  the  main  lobe  aperture  (angular  distance  between  main 
lobe  nulls)  is  equal  to  10°.  Determine  the  change  in  main  lobe  aperture 
due  to  a 10%  change  in  the  operating  frequency. 


12.  Consider  the  strip-line  of  figure,  of  known  dimensions  and  characteristic 

I M 

t§  X 


impedance  Zq.  The  strip-line  is  excited  with  dominant  TEM  mode  and  a pulsed 
voltage  V=V05(z-ct)  propagates  along  the  line.  Assume  the  electric  field 
distribution  to  be  uniform  along  x and  y in  the  strip- line  and  compute  the 
transient  far  field  in  the  E-plane  (x  = 0 plane) . Assume  the  thickness  of 
conductors  of  the  strip-line  to  be  infinitesimally  small. 
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